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ABSTRACT 

A  methodology  for  estimating  physical  parameters  in  a  class  of  structural  acoustic  sys¬ 
tems  is  presented.  The  general  model  under  consideration  consists  of  an  interior  cavity  which 
is  separated  from  an  exterior  noise  source  by  an  enclosing  elastic  structure.  Piezoceramic 
patches  are  bonded  to  or  embedded  in  the  structure;  these  can  be  used  both  as  actuators  and 
sensors  in  applications  ranging  from  the  control  of  interior  noise  levels  to  the  determination  of 
structural  flaws  through  nondestructive  evaluation  techniques.  The  presence  and  excitation 
of  the  patches,  however,  changes  the  geometry  and  material  properties  of  the  structure  as 
well  as  involves  unknown  patch  parameters,  thus  necessitating  the  development  of  parameter 
estimation  techniques  which  are  applicable  in  this  coupled  setting.  In  developing  a  framework 
for  approximation,  parameter  estimation  and  implementation,  strong  consideration  is  given 
to  the  fact  that  the  input  operator  is  unbonded  due  to  the  discrete  nature  of  the  patches. 
Moreover,  the  model  is  weakly  nonlinear  as  a  result  of  the  coupling  mechanism  between  the 
structural  vibrations  and  the  interior  acoustic  dynamics.  Within  this  context,  an  illustrating 
model  is  given,  well-posedness  and  approximation  results  are  discussed  and  an  applicable  pa¬ 
rameter  estimation  methodology  is  presented.  The  scheme  is  then  illustrated  through  several 
numerical  examples  with  simulations  modeling  a  variety  of  commonly  used  structural  acoustic 
techniques  for  system  excitation  and  data  collection. 
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1  Introduction 


The  recent  success  of  piezoceramic  materials  as  sensors  and  actuators  in  applications  involv¬ 
ing  structural  vibrations  has  spawned  intense  study  into  questions  regarding  the  modeling 
of  piezoceramic  actuator/sensor  interactions  with  underlying  structures,  the  optimal  design, 
placement  and  number  of  actuators/sensors  to  be  used,  and  the  development  of  effective  con¬ 
trol  strategies  in  a  variety  of  environments.  When  bonded  to  or  embedded  in  a  thin  structure 
(beam,  plate  or  shell),  piezoceramic  patches  derive  their  actuating  capacity  from  the  piez- 
electric  property  that  an  induced  voltage  produces  a  strain  in  the  material  thus  leading  to 
the  potential  for  producing  in-plane  forces  and/or  moments  when  the  patches  are  mounted  in 
pairs.  Conversely,  their  sensing  capabilities  are  due  to  the  dual  piezoelectric  effect;  namely,  a 
mechanical  force  leads  to  the  generation  of  a  proportional  voltage  across  the  element  which 
can  then  be  used  to  measure  accumulated  strain. 

However,  the  bonding  or  embedding  of  patches  in  the  underlying  structure  changes  not 
only  the  geometry  of  the  structure  but  also  physical  properties  such  as  the  density,  stiffness, 
Poisson  ratio  and  damping  coefficients,  when  indeed,  many  of  these  parameters  are  unknown 
even  for  the  homogeneous,  uniform  host  structure  material.  This  necessitates  the  development 
of  effective  parameter  identification  methods  to  be  used  when  estimating  system  parameters  in 
applications  such  as  experimental  model  validation,  the  determination  of  optimal  placement 
and  number  of  patches,  the  use  of  piezoceramic  patches  in  nondestructive  evaluation  (NDE) 
techniques,  and  the  implementation  of  model-based  control  schemes.  We  point  out  that  the 
estimation  of  physical  parameters  in  this  setting  differs  from  that  considered  in  much  of  the 
previous  literature  (see  [6]  and  the  references  therein)  in  that  here  the  patch  contributions 
to  the  system  lead  to  unbounded  (discontinuous)  input  and  output  operators  due  to  the  fact 
that  the  patches  cover  only  discrete  portions  of  the  structure. 

Parameter  estimation  methods  for  distributed  parameter  systems  involving  unbounded 
operators  have  been  developed  and  tested  in  the  case  in  which  piezoceramic  patches  are  used 
as  sensors  and  actuators  when  bonded  to  a  transversely  vibrating  beam  [15,  16].  There,  fit- 
to-data  techniques  involving  PDE  models  were  developed  which  could  be  used  to  estimate 
unknown  beam  parameters  given  various  data  forms.  Moreover,  in  that  setting,  results  per¬ 
taining  to  convergence  and  continuous  dependence  on  data  were  obtained  in  a  variety  of  cases 
involving  physically  tractable  methods  for  exciting  the  system  and  measuring  data. 

In  this  work,  we  develop  an  analogous  methodology  which  can  be  used  for  estimating 
physical  parameters  in  structural  acoustic  systems.  In  the  systems  of  interest,  an  exterior 
noise  source  is  separated  from  an  interior  cavity  by  a  thin  elastic  structure  (a  beam,  plate 
or  shell).  As  energy  is  transferred  from  the  exterior  field  to  the  structure,  vibrations  develop 
which  then  lead  to  unwanted  interior  noise  through  acoustic/structure  interactions.  Control  of 
this  unwanted  interior  noise  is  accomplished  through  sensing  and  actuating  via  piezoceramic 
patches  which  are  bonded  to  the  structure.  Before  model-bcised  control  schemes  can  be  imple¬ 
mented,  however,  the  physical  parameters  of  the  structure  (which  now  includes  the  patches) 
must  be  estimated  from  data  which  is  collected  both  on  the  structure  and  from  the  acoustic 
response  in  the  enclosed  cavity.  Although  similarities  exist  between  the  problem  of  estimat¬ 
ing  physical  parameters  for  the  isolated  structure  and  that  involving  the  structural  acoustic 
system,  the  hyperbolic  contributions  due  to  the  acoustic  component  sufficiently  change  the 


1 


problem  dynamics  so  as  to  warrant  in-depth  stud„  of  techniques  for  the  latter  coupled  system. 
Finally,  we  note  that  although  the  initial  impetus  for  developing  distributed  parameter  esti¬ 
mation  techniques  for  structural  acoustic  systems  was  motivated  by  model-based  noise  control 
considerations,  the  same  techniques  can  be  used  when  performing  vibration  analysis  or  using 
NDE  methods  to  determine  structural  flaws  in  these  coupled  systems. 

The  structural  acoustic  problem  used  here  to  motivate  and  illustrate  the  development  of 
an  appropriate  parameter  estimation  methodology  consists  of  a  2-D  enclosed  cavity  which  is 
separated  from  the  perturbing  exterior  noise  field  by  a  thin  beam.  This  model  represents 
a  2-D  slice  from  a  3-D  model  for  several  experimental  apparata  being  used  in  the  Acoustics 
Division,  NASA  Langley  Research  Center,  to  test  modeling,  parameter  estimation  and  control 
strategies.  We  add,  however,  that  the  methodolo'gy  being  presented  is  equally  valid  for  es¬ 
timating  parameters  in  many  3-D  models  representing  various  experimental  setups  currently 
in  use  (see  the  models  in  [9,  10]).  This  2-D  model  was  chosen  simply  because  it  simplifies 
the  discussion  and  more  clearly  illustrates  the  process  involved  in  developing  the  parameter 
identification  techniques. 

The  mode’  being  used  to  illustrate  the  methodology  is  weakly  nonlinear  due  to  the  manner 
in  which  the  structural  vibrations  couple  with  the  interior  acoustic  fields.  While  linearization 
provides  a  very  good  approximation  to  the  system  dynamics  (see  [1]),  we  retain  the  non¬ 
linearity  here  so  as  to  illustrate  some  of  the  general  analytic  assumptions  which  are  made 
when  extending  well-posedness  and  parameter  convergence  results  for  the  corresponding  lin¬ 
ear  problem  to  a  weakly  nonlinear  problem  of  this  type.  This  also  facilitates  the  demonstration 
of  numerical  techniques  which  can  be  used  when  simulating,  testing  and  implementing  the 
parameter  estimation  scheme  in  the  nonlinear  problem.  In  discussing  parameter  estimation 
methods  for  structural  acoustics  problems  of  this  type,  our  emphasis  is  on  the  formulation  of 
the  problem  in  a  manner  which  is  conducive  to  approximation  and  implementation  both  in 
the  linear  and  nonlinear  forms  cis  well  as  under  a  variety  of  damping  assumptions. 

In  the  second  section  of  this  presentation,  a  model  for  the  2-D  system  being  used  to 
illustrate  the  parameter  estimation  method  is  presented.  Details  regarding  the  modeling  of 
the  ax:oustic  and  structural  components  as  well  as  coupling  conditions  are  given,  and  care  is 
taken  to  motivate  the  assumptions  which  lead  to  various  damping  conditions  in  tlie  system 
model.  In  formulating  the  strong  form  of  the  system  model,  details  are  also  given  regarding 
the  interactions  between  the  piezoceramic  patches  and  the  underlying  structure  (beam)  as 
well  as  the  weakly  nonlinear  interactions  between  the  beam  and  the  interior  acoustic  field.  To 
provide  a  formulation  which  is  conducive  to  approximation  in  the  context  of  unbounded  input 
operators  as  well  cis  facilitates  the  discussion  of  well-posedness  results,  the  weak  form  of  the 
system  equations  is  then  developed  and  posed  in  terms  of  sesquilinear  forms  and  the  bounded 
operators  which  they  define.  Finally,  we  show  that  the  solution  trajectories  can  be  expressed 
in  terms  of  a  semigroup  on  an  appropriate  space  and  within  this  framework,  the  assumptions 
underlying  the  well-posedness  results  for  the  linear  and  nonlinear  problems  are  discussed. 

A  discretization  method  suitable  for  simulations  and  the  implementation  of  the  parameter 
estimation  method  is  outlined  in  Section  3.  This  discussion  is  kept  relatively  brief  since  details 
regarding  the  corresponding  finite  dimensional  system  for  the  linearized  problem  can  be  found 
in  [1].  However,  the  section  does  provide  additional  details  concerning  the  discretization  of  the 
nonlinear  component  of  the  operator  and  a  brief  algorithm  for  carrying  out  this  discretizetion 
is  included. 


2 


A  parameter  estimation  scheme  suitable  for  data  consisting  of  displacement,  velocity  or 
acceleration  mecisurements  on  the  beam,  voltage  measurements  from  the  patches,  or  pressure 
measurement  inside  the  cavity  is  presented  in  Section  4.  Assumptions  on  the  form  of  the 
unknown  parameters  are  discussed  and  conditions  leading  to  convergence  and  continuous 
dependence  on  data  results  for  the  linear  and  nonlinear  problems  are  outlined. 

In  Section  5,  numerical  examples  demonstrating  the  parameter  estimation  techniques  with 
a  variety  of  data  types  and  methods  for  exciting  the  system  are  presented.  Specifically, 
examples  are  given  in  which  the  force  to  the  system  is  provided  by  a  numerically  simulated 
acoustic  source,  simulated  voltage  inputs  to  the  patches,  and  a  simulated  hammer  impact 
to  the  beam.  The  simulated  data  under  comparison  in  the  examples  consists  of  acceleration 
values  of  the  beam,  voltage  (accumulated  strain)  values  from  the  patches,  and  voltage  values 
in  conjunction  with  interior  acoustic  pressure  values.  The  conclusions  from  this  study  and 
physical  considerations  concerning  the  implementation  of  the  method  are  summarized  in  the 
final  section  of  the  paper. 


2  System  Model 

The  model  of  interest  consists  of  an  exterior  noise  source  which  is  separated  from  an  interior 
cavity  U{t)  by  a  common  elastic  boundary  ro(<)  that  is  modeled  by  an  Euler- Bernoulli  beam 
as  depicted  in  Figure  1.  The  beam  is  assumed  to  have  length  a,  width  b  and  thickness  h. 
The  Young’s  modulus,  mass  density  (in  mass  per  unit  volume)  and  damping  coeflBcient  for 
the  homogeneous  beam  are  denoted  by  Eb,pb  and  cdj,  respectively.  Due  to  the  nature  of  the 
exterior  forces  and  the  manner  in  which  the  patches  are  excited,  we  will  be  considering  only 
transverse  vibrations  w{t,x). 

Bonded  to  the  beam  are  piezoceramic  patches  which  are  mounted  in  pairs  as  depicted 
in  Figure  1.  In  this  discussion,  it  is  assumed  that  the  patches  have  thickness  T,  width  6, 
Voung’s  modulus  Epe,  density  ppe,  and  damping  coeflScient  cpp^.  Moreover,  it  is  assumed 
that  the  bonding  layers  for  each  patch  have  equal  thicknesses.  Young’s  moduli,  densities  and 
damping  coefficients,  and  these  parameters  are  denoted  by  T^,  pbe  and  cubi-,  respectively. 
We  emphasize  that  these  assumptions  are  made  solely  for  clarity  of  presentation,  and  similar 
results  can  be  obtained  in  an  analogous  manner  for  the  more  general  case  in  which  the  patches 
and  bonding  layers  have  differing  thicknesses  and  material  properties  (see,  for  example,  [12]). 

Finally,  we  assume  that  inside  the  cavity,  there  is  a  region  =  IJill  which  provides 
a  first  approximation  to  the  interior  objects  (eg.,  passengers,  seats,  storage  compartments) 
which  disrupt  and  damp  the  interior  acoustic  fields.  This  region  is  assumed  to  have  positive 
measure  and  is  taken  to  be  small  as  compared  to  n(t)  (see  Figure  1  where  Nr  =  4). 

2.1  Acoustic  Component 

We  consider  first  the  acoustic  wave  dynamics  inside  the  cavity  Q,{t).  The  variables  of  interest 
to  us  are  the  pressure  P,  density  p,  and  the  velocity  U,  each  of  which  can  be  represented  in 
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Fig.  1.  The  2-D  domain. 


terms  of  a  mean  and  fluctuating  component 

P{t,x)  =  Po{x)  +  p{t,x) 

pityx)  =  po{x)  +  p{t,x) 

(j{t^x)  =  Oo{x)  +  u{t,x) 

(here  we  have  taken  x  =  (x,  y)  and  are  assuming  that  the  rate  of  sound  travel  is  sufficiently 
rapid  so  that  little  heat  transfer  takes  place).  For  the  range  of  magnitudes  involved  in  these 
problems  (<  140  dB),  it  is  usual  to  eissume  linear  relations  when  considering  constitutive  laws 
and  force  balancing  [20],  and  we  make  that  assumption  throughout  the  analysis  which  follows. 

In  the  region  air  damping  is  omitted  due  to  the  relatively  small  dimensions  of  the 

type  of  experimental  cavities  being  simulated.  Hence  in  that  region,  an  increase  in  pressure 
brings  about  a  proportional  strain  with  the  ratio  defined  as  the  bulk  modulus  of  elasticity  Ej 
where 

=  or  p=-E,V-s  (2.1) 

{s{t^x)  denotes  the  displacement  of  the  center  of  gravity  of  an  infinitesimal  element  of  the 
medium  and  satisfies  St  =  u).  On  the  other  hand,  the  maierial  objects  lumped  in  Cl  will 
provide  sorhe  medium  damping  and  here  we  assume  that  a  change  in  pressure  yields 

p  =  -EV  •  s  -  dV  •  i;  (2.2) 

where  E  and  d  denote  the  bulk  modulus  of  elasticity  and  damping  coefficient  of  the  medium 
in  0.  We  point  out  that  this  use  of  a  generalized  Hooke’s  law  in  which  stress  is  proportional  to 


a  linear  combination  of  strain  and  strain  rate  is  done  under  the  assumption  that  relatively  low 
acoustic  frequencies  are  excited  (<  1000  Hz),  and  is  similar  to  the  constitutive  law  leading  to 
Kelvin- Voigt  damping  of  vibrations  in  elastic  materials.  We  also  emphasize  that  this  damping 
model  should  be  considered  as  a  first  approximation  to  the  actual  acoustic  damping  mechanism 
in  the  medium  Ct,  and  depending  on  the  specific  materials  involved,  the  manner  of  acoustic 
excitation,  and  the  geometry  of  the  physical  system,  more  comprehensive  models  may  be 
required  to  accurately  describe  the  medium  damping. 

Force  balancing  in  the  acoustic  cavity  yields  the  relation  (the  Euler  equation) 

Po^  =  -Vp  (2.3) 


where 


(  pf  ,  X  e  f)(t)/n 

\  p  ,  £  e  fi 


denotes  the  equilibrium  density  of  the  medium  ((2.3)  is  equivalent  to  the  linearized  momentum 
equation  in  fluid  dynamics).  By  taking  two  time  derivatives  of  (2.1)  and  (2.2),  the  divergence 
of  (2.3),  and  eliminating  cross  terms,  one  arrives  at 


dP  PS  ^ 

—  =  -Ap-b-zApt 

ar  p  p 


X  € 


£  €  fi . 


Taking  the  curl  of  the  mome..tum  equation  (2.3)  and  noting  that  u  =  St  we  obtain  |(Vx 
Pou)  =  0.  Hence  the  vorticity  w  =  V  x  pou  is  constant  in  time  in  Sl{t).  Under  the  assumption 
that  the  initial  vorticity  w(0)  is  zero,  we  may  conclude  that  V  x  pou  =  0  for  all  time  or  that 
the  flow  is  irrotational  in  fi(t).  Thus  in  fl(t),  there  exists  a  scalar  velocity  potential  ^  such 
that 

pou  =  -V<l>.  (2.5) 

The  flow  is  more  complex  near  and  in  the  region  as  a  result  of  the  viscous  effects  and 
medium  damping.  This  can  potentially  lead  to  rotational  components  in  the  acoustic  field 
which  in  general  necessitates  the  use  of  a  vector  potential.  As  a  first  approximation,  however, 
we  are  assuming  that  the  rotational  components  near  and  in  0  are  negligible  and  a  relation 
of  the  form  (2.5)  will  be  used  throughout  the  acoustic  cavity  fi(t). 

Substituting  (2.5)  into  the  momentum  equation  (2.3)  we  find 


which  implies  that  p  =  since  no  acoustic  sources  or  sinks  are  present  in  fl(t).  Use  of  this 

at 

pressure  expression  (actually  pt  =  (f>tt)  after  differentiation  in  time  once  in  equations  (2.1), 


(2.2)  then  yields 


dt^  PS  ^ 

-7—  =  —A6+-A6t 


X  6  0,{t)/Cl 


xen . 


Finally,  liardwall  boundary  conditions  are  assumed  since  the  experimental  apparatus  being 
modeled  heis  concrete  walls.  With  n  denoting  the  outward  unit  normal  to  the  cavity  and  c^(x) 
and  d{x)  denoting  the  speed  of  sound  and  damping  coefficients  g.  ■ni  by 


c^(x)  —  /  ®  ^  Cl{t)/^ 

1  Ejp  ,  X  €  n 


d{x)  = 


0 

djp 


X  €  n{t)/Q 

X  €  , 


the  acoustic  response  can  be  modeled  cis 


dt^ 


=  c^A<l>  +  dA(l>t 


{x,y)  6  n(t)  ,^  >  0  , 


-11  =  0 


(x,i/)  G  r  >  0 . 


(2.6) 


We  emphasize  that  this  model  was  derived  under  the  assumption  that  the  only  acoustic 
damping  in  the  cavity  occurs  in  the  region  0  and  hence  d{x)  =  0  in  the  rest  of  the  chamber, 
i.e.  for  X  G  Moreover,  we  have  assumed  that  the  flow  is  irrotational  in  the  region 


2.2  Beam  Component 

Through  force  and  moment  balancing,  the  equation  of  motion  for  the  transverse  displacements 
w  of  the  beam  are  found  to  be 


d^M 

pwtt  +  -^{t,x)  =  -^t(t,x,u;(t,x))H-/(<,x) 

w(t,0)  =  0 


x0<x<a,  t>0, 
i>0  , 


where  M  is  the  total  beam  moment,  /  is  the  force  due  to  pressure  from  the  exterior  noise  field 
and  4t{t,  x,  w{t,  x))  is  the  backpressure  due  to  the  ensuing  acoustic  waves  inside  the  cavity 
(this  latter  term  is  in  general  nonlinear  since  its  effect  occurs  on  the  surface  of  the  vibrating 
beam).  For  pairs  of  patches  having  edges  xi  and  X2,  the  density  of  the  structure  is 


p{x)  =  pbhb  +  2b{pbeTu  +  ppeT)  Xpeix) 


(2.7) 


where  the  characteristic  function 

Xpe(a:)  = 


1 

0 


Xi  <  X  <  X2 
otherwise 


localizes  the  patch  effects  between  the  endpoints  xi  and  X2  (see  [11,  12]  for  details). 
The  general  beam  moment 


M{t,x)  =  M{t,x)  —  Mpc{t,x) 

consists  of  an  internal  component  Af,  depending  on  material  and  geometric  properties  of  the 
beam  and  patches,  and  an  external  component  Mpe  (the  control  term)  which  results  from  the 
activation  of  the  patches  through  an  applied  voltage  (see  Figure  2). 
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For  a  beam  undergoing  pure  bending  motion  with  out-of-phase  excitation  of  the  patches, 
the  inteinal  and  external  moments  are  given  by 

XML  \  ,-7/  77  \ 

A/(i,x)  =  £/(.)^  +  cc/(x)— 

My^{t,x)=K^VXr,e{x) 

where  V  is  the  voll-*je  into  the  patches.  As  shown  in  [11],  the  stiffness,  damping  and  control 
constants  h  -r  the  combined  structure  are 

h^h  26  r  ,  , 

EI[x)  =  Eb—  +  y  [Ebta^bi  +  ipeOSpe]  Xpe(®) 

CdI (x)  =  C£)6^  +  y  [cDhld-m  +  CDpeOSpe]  Xpe(x) 

—  Epebd^i^h  -f  2Tbe  +  T)  . 

The  ci/D-^tants  036^  and  given  by  asbt  =  {hf2  +  Tbef  -  {h/2y  and  oape  =  {h/2  +  Tbe  +  T)^- 
(h/2  -f  Ti  )  result  from  the  integration  of  stresses  through  the  bonding  layer  and  patch,  and 
dji  a  pi'  uceramic  constant  which  relates  the  amount  of  strain  produced  in  the  patch  to 
the  h'v't'i  yj'i.  voltage  being  applied. 

no  tais  point  it  is  worth  commenting  further  on  the  damping  term  cope  which  is  taken  to 
be  a  combination  of  the  Kelvin  Voigt  damping  coefficient  for  the  patch  and  the  damping  which 
results  from  the  production  of  current  when  the  structure  vibrates.  This  latter  contribution  to 
the  damping  results  from  the  piezoelectric  effect  of  the  patches  which  dictates  that  a  voltage 
is  pr(/'luced  when  the  patch  is  subjected  to  in-plane  strains.  Under  the  eissumption  that  the 
Kelvin- Voigt  (material)  and  electrical  damping  have  approximately  the  same  types  of  effect 
in  the  patch,  we  have  combined  the  two  into  the  coefficient  copg  which  must  be  considered  to 
be  unknown  and  like  the  other  parameters,  must  ultimately  be  estimated  using  data  fitting 
techniques  with  experimental  data  when  considering  actual  applications.  We  also  point  out 
that  the  expression  (2.8)  can  easily  be  generalized  to  include  the  possibility  of  differing  material 
properties  in  the  two  patches  or  bonding  layers  (again,  see  [11,  12]). 

The  fact  that  the  patches  generate  a  voltage  when  strained  implies  that  they  can  be  used 
as  sensors  Jis  well  as  actuators.  Specifically,  the  voltage  produced  by  the  patches  when  the 
beam  bends  is 


V,{t)  =  IC^  f 

Jx 


d^w , 


(/.,  x)dx 


dw .  ,  dw , 


where  the  sensor  constant  depends  upon  piezoelectric  material  properties  as  well  as  the 
geometi  '  and  size  of  the  patch  (see  [18]).  Hence  the  voltage  provides  a  measure  of  accumulated 
strain  in  the  beam  (see  Figure  2)  thus  enabling  the  patch  to  serve  as  a  sensor  in  a  variety  of 
applications  involving  the  measurement  of  beam  vibrations. 
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Fig.  2.  Acoustic  cavity  with  piezoceramic  patches;  (a)  fixed  domain  ft  and  bouridary  F  U  Fo;  (b)  patch 
generation  of  pure  bending  moments,  and  (c)  patch  sensing  of  accumulated  strain. 

2.3  Coupling  Conditions 

In  the  model  discussed  thus  far,  the  structural  and  internal  acoustic  responses  are  coupled 
thiough  the  backpressure  (f>t{t,x,w{t,x))  on  the  surface  of  the  beam.  A  second  coupling 
equation  is  the  continuity  of  velocity  (or  momentum)  condition 

wAt,  x)  =  —  V<^(i,  X,  w{t,  x))-  j  ,  0<x<a,t>0 

PJ 

which  results  from  the  assumption  that  the  beam  is  impenetrable  to  air.  We  point  out  that  the 
velocity  condition  provides  a  form  of  damping  to  the  beam  which  is  similar  to  that  obtained 
with  the  incorporation  of  viscous  (air)  damping  elfects  (modeled  by  a  term  of  the  form  'ywt  in 
the  beam  equation).  As  noted  in  the  examples,  the  internal  Kelvin- Voigt  damping  in  concert 
with  the  coupling  effects  due  to  the  continuity  of  velocity  and  backpressure  causes  a  beam 
response  which  differs  somewhat  from  that  observed  with  an  uncoupled,  undamped  beam 
having  the  same  dimensions. 

The  model  we  have  developed  has  nonlinearities  in  the  (i)  variable  domain  ft(i),  (ii)  back 
pressure  term  —(j)t{t,x,w{t,x)),  and  (iii)  velocity  term  ■^'^<f>{t,x,w{t,x))  ■  j.  Under  an  as¬ 
sumption  of  small  displacements  {w{t,x)  =  w{t,x)  -f  6  where  u)  =  0)  which  is  inherent  in 
the  Euler- Bernoulli  formulation,  the  variable  domain  fl{t)  can  be  approximated  by  the  fixed 
domain  0  =  [0,o]  x  [0,£]  as  shown  in  Figure  2.  Note  t^'at  with  this  cissumption,  the  velocity 
term  —'V(j){t,x,w{t,x))  •  j  can  be  approximated  by  the  normal  term  —’V(l){t,x,w{t,x))  •  ii 
which  arises  when  developing  the  weak  form  of  the  equation.  The  fully  nonlinear  form  of  the 
back  pressure  coupling  term  is  retained  throughout  the  following  discussion. 


2.4  Strong  Form  of  the  System  Equations 

For  the  coupled  system  in  which  s  pairs  of  patches  are  bonded  to  the  beam  and  excited  out- 
of-phase,  the  acoustic,  structural  and  coupling  components  just  discussed  can  be  combined  to 
yield  the  approximate  nonlinear  model 

=  c^A(f)  -1-  dA^t  {x,y)  6  Q,  ,t  >  0  , 

S/<j)  •  n  =  0  (x,  j/)  6  r  ,  f  >  0  , 

X,  w{t,  x))  •  h  =  wAt,  x)  0<x<a,<>0, 

PS 


f  d^w  d^w 


=  X,  w{t,  x))  -1-  /(t,  x)  +  —j2  }^fui{i)Xve.{x) 

w{t,0)  = =  w{t,a)  = =  0  t>0  , 

^(0,  X,  y)  =  (l)o{x,  y)  ,  tt;(0,  x)  =  wo{x) 

3!,  y)  =  (f>x (x,  y)  ,  u;e(0,  x)  =  lUi (x)  . 


0  <  X  <  a  , 
t  >  0  , 


(2.10) 


Here  Ui{t)  is  the  voltage  being  applied  to  the  patch  and  Xpe,  denotes  the  characteristic 
function  over  the  patch. 

We  point  out  that  the  piezoceramic  material  parameters  Kf  (see  (2.9))  and  /Cf ,  i  = 
as  well  as  the  beam  parameters  p,EI  and  cpI  are  considered  to  be  unknown  <nd 
are  estimated  using  inverse  problem  techniques  ais  discussed  in  later  sections.  While  the 
expressions  given  in  (2.7)  and  (2.8)  can  be  used. as  starting  values  in  the  parameter  estimation 
routines,  experimental  evidence  (see  [15,  16])  has  indicated  that  the  final  parameter  values 
can  vary  quite  significantly  from  the  analytic  values  due  to  the  contributions  from  the  bonding 
layer,  variation  in  the  measurement  of  physical  constants,  and  nonuniformities  in  the  various 
materials.  This  combined  with  the  lack  of  analytic  expressions  for  the  damping  constant 
necessitates  the  estimation  of  these  parameters  before  model-based  control  strategies  can  be 
implemented. 

We  also  emphasize  that  the  parameters  K.f,p,EI  and  cpI  are  piecewise  constant  in  na¬ 
ture  due  to  the  presence  and  differing  material  properties  of  the  bonding  layer  and  patches 
(see  (2.8)  as  well  as  the  results  in  [15]).  This  leads  to  difficulties  with  the  strong  form  of  the 
system  equations  since  it  necessitates  the  second  derivatives  of  the  Heaviside  function  (equiv¬ 
alently,  derivatives  of  the  Dirac  delta)  thus  yielding  an  unbounded  control  input  operator. 
The  differentiation  of  the  discontinuous  material  parameters  also  leads  to  difficulties  when  ap¬ 
proximating  the  dynamics  of  the  coupled  system.  To  avoid  these  problems,  it  is  advantageous 
to  formulate  the  problem  in  weak  or  variatic.’ial  form  (the  use  of  the  variational  form  also 
permits  the  use  of  beisis  functions  having  less  smoothness  than  required  for  those  used  when 
approximating  the  solution  to  the  strong  form  of  the  equations). 


Finally,  we  note  that  in  the  Ccise  of  no  acoustic  cavity  damping  {d  =  0),  the  model  (2.1G) 
is  completely  equivalent  to  the  nonlinear  models  that  are  the  basis  of  the  investigations  in 
[1,  8,  10,  11]  if  one  replaces  the  coupling  terms  pj(l)t  and  v  =  —  V<^  in  those  models  by  and 
V  =  — respectively.  That  is,  the  potential  used  in  those  references  differs  from  the  one 
used  here  by  a  multiplicative  factor  — . 


2.5  Weak  Form  of  the  System  Equations 

An  appropriate  choice  for  the  state  of  the  second-order  problem  (2.10)  is  the  pair  (<^, 
consisting  of  the  acoustic  potential  and  beam  displacement.  It  follows  immediately  that  with 
this  choice,  the  state  for  the  problem  in  first-order  form  \s  Z  =  (<^,  tu,  w)  which  contains  the 
pressure  (since  p  =  ^)  as  well  cis  the  beam  displacement  and  velocity. 

The  state  space  and  space  of  test  functions  are  taken  to  be  the  product  spaces  7i  =  V  x  H 
and  V  =V  xV  where  the  Hilbert  spaces  H  and  V  are  given  hy  H  =  Z-'^(fl)  x  L'^(ro)  and  V  = 
X  ffg(ro).  Here  Z^(Q)  and  H^(Q,)  denote  the  quotient  spaces  of  and  over  the  con¬ 
stant  functions  and  Hq^Tq)  is  given  by  ^^0(^0)  =  {0  €  //^(Fo)  :  fpix)  =  ip'{x)  =  0  at  x  =  0,a} 
(the  use  of  the  quotient  space  results  from  the  fact  that  the  potentials  are  determined  only 
up  to  a  constant).  From  energy  considerations,  the  H  and  V  inner  products  are  taken  to  be 


((»)’(  9  ^  /r. 


while  the  product  space  inner  products  are  given  by 


((*),()[))  =(*.T)v +  ('!', A)^  ■ 


Integration  by  parts  in  a  manner  analogous  to  that  in  [1]  in  combination  with  the  approx- 
iiriation  of  the  variable  domain  n(t)  by  the  fixed  domain  =  [0,  a]  x  [0,  and  the  use  of 
Green’s  theorem  then  yields  the  nonlinear  first-order  variational  form 


+  f  /  El D^wD^pd-y 

»/n  Jtq 

+J^  -^d  ^cdI +  [^(lo)?/  —  |  d-y 

=  +  f'nd'y 
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for  all  ^  =  (^,  rj)  G  V.  We  point  out  that  in  this  variational  form  the  derivatives  have  been 
transferred  from  the  plate  and  patch  moments  onto  the  test  functions.  This  eliminates  the 
problem  of  having  to  approximate  the  derivatives  of  the  characteristic  function  and  the  Dirac 
delta  eis  is  the  case  with  the  strong  form  of  the  equations. 

The  system  (2.11)  can  be  formally  approximated  by  replacing  the  state  variables  by  their 
finite  dimensional  approximations  and  constructing  the  resulting  matrix  system,  and  it  is  in 
this  form  that  we  will  consider  approximation  strategies  in  Section  4.  In  order  to  discuss  the 
well-poseduess  of  the  problem,  however,  it  is  advantageous  to  pose  the  problem  in  an  abstract 
Cauchy  format  eis  discussed  in  the  next  subsection. 

2.6  Abstract  First  Order  Formulation 

Following  the  theoretical  work  in  [2,  3],  it  is  advantageous  to  formulate  the  problem  in  terms 
of  sesquilinear  forms  and  the  bounded  operators  which  they  define  (see  also  [1,  10]  for  further 
examples  detailing  the  abstract  formulation  of  structural  acoustic  systems  in  this  manner). 
We  begin  by  pointing  out  that  the  Hilbert  spaces  H  and  V  form  a  Gelfand  triple  K  ~ 

H*  V*  with  pivot  space  H  (further  details  concerning  the  basic  definitions  and  fundamental 

functional  analysis  theory  here  can  be  found  in  [22]).  For  =  (^,  lo)  and  ^  =  (^,??)  in  V,  we 
then  define  the  sesquilinear  forms  a-,  :V  xV—*C,  i  =  1,2  by 

ai($,  <p)  =  /  V«^  •  Verfw  +  f  EID^wD^rid'y  , 

^  (2.12) 
cr2($, «»)  =  /  +  /  {cdID^wD^tj  +  -  t«^]}  d'y  . 

With  these  definitions,  it  is  straightforward  to  show  that  and  (T2  satisfy  various  conti¬ 
nuity,  symmetry  and  coercivity  conditions.  Namely,  cti  satisfies 

|cri($, '5)1  <  ci|4»|v|'5|v'  ,  for  some  ci  6  IR  (bounded) 

Re  (ri(^>,$)  >  02]$]^  ,  for  some  C2  >  0  (l/-elliptic)  (2.13) 

(Ji($,  ^)  =  ai('5,  ^►)  (symmetric) 

for  all  €  V  f  .he  boundedness  results  from  Schwarz’s  inequality  for  inner  products  in 
conjunction  with  equivalence  results  for  various  Sobolev  norms,  while  the  K-ellipticity  and 
symmetry  of  cri  follow  directly  from  the  fact  that  <ri($,  ^)  =  ($,  ^)y). 

The  behaviour  of  (T2  depends  upon  the  form  of  damping  in  the  cavity  0.  If  damping  is 
included  in  a  region  ft  having  positive  measure  as  discussed  in  Section  2.1,  then  cr2  satisfies 


|cr2($,^)|  <  cal^'lvj'^jv  ,  for  some  C3  €  IR  ^bounded) 
Re  (72($,  $)  >  C4l4>|y  ,  for  some  C4  >  0  (V-elliptic)  . 


(2.14) 


The  V-ellipticity  follows  from  the  fact  that 


Re  (r2($,  $)  =  [  \d  \Vcf>fduj  +  f  cdI  {d^wV  d'y 

V  fi  c  v  r^o 

>  h  \V(l>\^diJ  +  k2j^  El  [d'^w)  '  (i7 

>  cM 

with  the  final  inequality  resulting  from  the  observation  that  there  exists  a  constant  such 
that  ki  /q  \V(f>^dit}  >  ks  |V<^P<L;  for  all  4>  €  V.  On  the  other  hand,  if  the  cavity  is  taken 
to  be  completely  open  with  no  medium  damping  assumed  (hence  mea.s(l))  =  0  and  d  =  0 
throughout  Q),  then  one  can  only  establish  the  weaker  //-semiellipticity  condition 

Re  a-2($,  >  C4|$|^  ,  C4  >  0 

for  <r2  (in  this  case  C4  =  0  since  (T2($,4>)  contains  no  acoustic  components  when  damping  is 
omitted  in  the  cavity).  As  discussed  in  the  next  section,  the  inclusion  of  medium  damping 
within  the  subregions  Q  and  the  resulting  K-elliptic  behavior  of  <72  yields  stronger  semigroup 
results  than  those  obtained  when  acoustic  damping  is  omitted. 

To  account  for  the  patch  contributions  when  a  voltage  is  applied,  we  let  U  denote  the 
Hilbert  space  containing  ^he  voltage  inputs  and  we  define  the  operator  B  €  C{U,  V*)  by 

(Bu,^)y.  v=  j  ^/Cfu,Xpe,-£>V7 

for  €  K,  where  {’,')vy  is  the  usual  duality  pairing.  Finally,  the  external  forcing  and 
nonlinear  perturbation  terms  are  given  by  F  =  {0,flpb)  and  G{z,zt)  =  (0,  — ^((iw))  where 
^^t(w)  =  ^i{t,x,w{t,x))  =  (j)t{t,x^w{t,x))  —  ^t(t,a;,0)  denotes  the  nonlinear  perturbation  to 
the  linear  coupling  term. 

With  these  definitions,  we  can  write  the  system  (2.11)  in  the  abstract  weak  or  variational 
form' 


Mt),  ^)v,v  +  ^)  +  =  {Bu{t)  +  F  +  G{z,  zt),  'I')v,.  V.  (2.ir3) 

for  in  V.  We  reiterate  that  the  state  for  the  problem  in  second-order  form  is  given  by 
z(t')  =  mV  H. 

To  pose  the  system  in  first-order  form,  we  form  the  product  space  terms  Bu{t)  =  (0,  Bu{t)), 
T{t)  —  (0, F(t))  and  Q{Z{t))  =  {Q^G{z{t),zt{t)))  in  V*  =-V  xV*  and  define  the  operators 
Au  V*)hy 

for  z  =  1,2  (the  existence  of  A\  and  A2  is  guaranteed  by  the  boundedness  of  cti  and  (72). 
Then,  for  the  state  Z{t)  =  {z{t),zi{t))  =  (^,to, ^,io)  m'H=VxH,  the  weak  form  (2.15)  is 
formally  equivalent  to  the  system 

Zt{t)  =  AZ{t)+C{t,Z{t))  (2.16) 
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in  V*  where 


C{t,  Z{t))  =  Bu{t)  +  r{t)  +  g{Z{t)) 


(2.17) 


and 

dom  ^  =  {0  =  (T,  A)  €  :  A  G  V,  /liT  +  AjA  €  H] 

r  0  /I  (2-18) 

A  = 

—A\  —A2  _ 

(further  details  concerning  the  formulation  of  the  first-order  system  in  the  linear  case  can 
be  found  in  [1]).  The  representation  (2.16)  is  formal  in  the  sense  that  the  manner  in  which 
differentiation  and  the  resulting  solution  exists  has  yet  to  be  specified.  This  will  be  discussed 
next  in  the  context  of  proving  well-posedness  results  for  the  system  model. 

Finally,  in  discussing  well-posedness  results  and  parameter  convergence  in  the  sections 
which  follow,  it  proves  useful  to  compare  the  nonlinear  system  with  that  which  results  when 
both  coupling  terms  are  linearized.  The  latter  is  found  by  replacing  the  term  ^t(t,  x,  w)  by 
<^((t,x,0)  in  the  various  expressions  for  the  coupled  system  (equivalently,  take  G{z^Zi)  =  0). 
This  yields  the  second-order  variational  form 

{zttii),  '^)v,V  +  0-2(2t(t),  'I')  =  {Bu{t)  -f  F,  (2.19) 

for  ^  €  V  and  consequently  the  first-order  system 

Zt{t)  =  AZ{t)-\-Bu{t)  +  J=^{t)  (2.20) 

in  V*. 

2.7  Model  Well-Posedness 

The  first  step  in  determining  the  well-posedness  of  the  semilinear  system  model  is  to  argue 
that  A  generates  a  semigroup  on  W.  As  noted  earlier,  the  sesquilinear  form  ctj  is  V-elliptic, 
continuous  and  symmetric  while  is  continuous  and  F-elliptic  if  damping  is  included  in  the 
region  Q  or  /f-semielliptic  if  damping  is  omitted  in  the  cavity.  In  both  cases,  the  Lumer- 
Philips  theorem  (with  arguments  similar  to  those  given  in  pages  82-84  of  [2])  can  be  used  to 
prove  that  the  operator  A  defined  in  (2.18)  generates  a  Co-semigroup  T  on  the  state  space  7i. 
The  semigroup  satisfies  the  exponential  bound  \T{t)\  <  for  t  >  0  (where  in  fact,  u  —  0  due 
to  the  fact  that  A  is  dissipative  as  shown  in  [2]).  Moreover,  if  medium  damping  is  included 
in  the  region  (which  implies  that  (T2  is  F-elliptic),  the  semigroup  T  is  analytic  on  H. 

For  the  problem  thus  posed,  the  state  lies  in  H  which  implies  that  the  semigroup  T 
generated  by  A  :  dom  A  CH-^  His  defined  on  H.  The  nonlinear  forcing  term  C{t,Z{t))  = 
Bu{t)-\-J^{t)  +  Q{Z{t)),  however,  lies  in  V*  rather  than  H  since  the  control  term  B  G  C{U,  V*) 
defines  the  product  space  control  term  Bu{t)  =  (0,  Bu(t))  €  {0}  x  K*  C  F  x  K*  =  V*. 
To  remedy  this,  “extrapolation  space”  ideas  and  arguments  similar  to  those  presented  in 
[3,  4,  19]  are  used  to  extend  the  semigroup  T{t)  on  to  a  semigroup  T{t)  on  a  larger  space 
Vy*  D  {0}  X  F*  so  eis  to  be  compatible  with  the  forcing  term. 


As  detailed  in  [10],  the  space  of  interest  is  defined  in  terms  of  dom  A*  where 
dom  A*  =  {x  =  (^>,  G  €  V,  AjtE*  -  6  H} 

(  -'I' 

A*y  =  I 

^  V 

Specifically,  the  space  W  =  [dom  A*]  is  taken  to  be  dom  A*  with  the  inner  product 

($,^)^  =  ((Ao-^*)<[>,(Ao-A*W^ 

for  some  arbitrary  but  fixed  Ao  with  Aq  >  a;  (recall  that  the  original  solution  semigrouj) 
satisfies  the  bound  |T(t)|  <  e*"*).  As  proven  in  [4],  the  resulting  W  norm  is  equivalent  to 
the  graph  norm  corresponding  to  A*.  Moreover,  we  have  that  {0}  x  V*  C  W*  =  [dom  A*]* 
(see  [10]  for  details). 

From  the  definition  of  A*  and  the  equivalence  of  the  VV  norm  with  the  graph  norm  corre¬ 
sponding  to  A*,  we  can  define  AQ  G  W*  by 

(i0)  (x)  = 

for  all  0  e  H,  X  G  W.  With  this  definition  and  the  Riesz  representation  theorem,  it  is 
shown  in  [10]  that  A  is  an  extension  of  the  original  operator  A  from  dom  A  CH  to  all  of  H. 
Finally,  as  proven  in  [4],  the  operator  A  is  the  infinitesimal  generator  of  a  Co-semigroup  '!’(<) 
on  W*  which  is  an  extension  of  T (t)  from  H  to  W*  (note  that  T (t)  is  also  analytic  if  medium 
damping  is  included  within  the  region  Cl). 

Having  extended  the  operator  A  and  hence  the  generated  semigroup  to  a  space  which  is 
compatible  with  the  forcing  function,  we  are  now  in  a  position  to  discuss  criteria  on  the  input 
terms  F  and  B  which  guarantee  the  existence  of  a  unique  solution  to  the  system  model.  In 
the  corresponding  linear  problem,  under  reasonable  regularity  conditions  on  t  i->  u{t)  and 
t  i-t  F{i),  one  can  immediately  argue  the  existence  of  a  unique  strong  solution  to  the  system 
in  terms  of  the  extended  semigroup  T{t).  For  the  semilinear  problem  of  interest,  however,  the 
nonlinear  nonhomogeneous  terms  must  satisfy  certain  continuity  criteria  in  order  to  obtain 
similar  results.  For  example,  if  we  let  X  denote  the  reflexive  Banach  space  W*  and  assume 
that  C  :  [0,r]  X  Af  — >  AT  defined  in  (2.17)  is  continuous  in  t  on  [0,  T]  and  uniformly  Lipschitz 
continuous  on  X,  then  the  integral  equation 

Z(t)  =  tws.  + 1  tit  -  s)  (  ^  ^  )  ..  (2.21) 

is  well-defined  for  Bu  +  F  +  G{Z)  G  L^((0,r),  V*).  Moreover,  for  Z(0)  =  Zq,  the  solution 
Z{t)  of  (2.21)  is  a  unique  mild  solution  to  (2.16)  (see  Theorem  1.2,  page  184  of  [21]).  In 
addition,  if  C  :  [0,  T]  x  X  — »  Af  is  Lipschitz  continuous  in  both  variables,  then  it  follows  from 
Theorem  1.6,  page  189  of  [21]  that  (2.21)  provides  the  strong  solution  to  (2.16)  interpreted  in 
the  Vy*  sense. 

The  required  continuity  of  the  nonhomogeneous  terms  Bu  and  F  is  demonstrated  in  [10] 
and  hence  the  remaining  question  concerns  the  Lipschitz  continuity  of  the  nonlinear  coupling 
term  G{z,zt)  =  (0,  — <^t(t«))-  If  we  assume  that  the  input  terms  F  and  Bu  are  sufficiently 
smooth  so  as  to  assure  the  necessary  continuity  in  G{z,zt),  then  the  nonlinear  system  is 
well-posed.  These  results  are  summarized  in  the  following  result. 
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Lemma  1.  (Well-Posedness  of  the  Nonlinear  System)  Consider  the  nonlinear  system 
represented  by  (2.15),  (2.16)  or  (2.21)  and  assume  that  F  and  Bu  are  sufficiently  smooth 
so  that  C{i,Z{t))  =  Bu{t)  +  F{t)  +  Q{Z{t))  is  Lipschitz  continuous  in  both  t  and  Z. 
Then  (2.21)  with  Zq  E  H  \s  the  unique  strong  solution  to  (2.15)  both  when  acoustic 
damping  is  eissumed  in  the  region  Cl  and  when  no  damping  is  present  in  the  cavity  (in 
which  case,  <72  is  only  //-semielliptic). 

These  results  can  be  further  extended  when  acoustic  damping  is  assumed  in  the  region  Cl 
and  hence  (Ja  is  K-elliptic  and  T  is  analytic.  In  this  case,  the  mild  (and  thus  strong)  solutions 
to  (2.20)  are  guaranteed  to  be  equivalent  to  the  weak  or  variational  solutions  to  (2.19)  for 
sufficiently  smooth  nonhomogeneous  terms.  The  following  theorem  summarizes  this  result  for 
the  linearized  problem. 

Theorem  1.  (Linearized  System:  Equivalence  of  Solutions)  Consider  the  system  rep¬ 
resented  by  (2.19)  or  (2.20)  and  suppose  that  the  mappings  t  u{t)  and  1 F{t)  from 
[0,  T]  to  IR^  and  V*,  respectively,  are  Lipschitz  continuous.  Furthermore,  assume  that 
medium  damping  of  the  form  (2.2)  is  present  in  the  region  Cl  C  0,.  Then  for  each 
Zo  €  H  =  dom  A,  we  have  that  (2.20)  taken  with  Z{0)  =  Zq  has  a  unique  strong 
solution  given  by  (2.21)  with  G{2)  =  0.  Moreover,  this  strong  solution  is  equal  to  the 
weak  solution  of  (2.19). 

The  proof  of  the  equivalency  between  strong  and  weak  solutions  follows  that  given  in  [3] 
for  general  second-order  systems  with  unbounded  input  terms.  We  point  out  that  numerous 
numerical  results  have  indicated  similar  results  for  the  nonlinear  problem  and  the  case  in  which 
damping  is  omitted  in  the  cavity  even  though  we  have  not  extended  the  results  equating  the 
strong  and  weak  solutions  to  cover  those  cases. 


3  State  Approximation 

The  modeling  and  well-posedness  discussion  thus  far  has  been  for  the  infinite  dimensional  non¬ 
linear  structural  acoustic  system.  In  this  section  we  discuss  a  Galerkin  scheme  for  discretizing 
the  problem  which  can  be  used  when  simulating  the  system  dynamics,  estimating  the  physical 
parameters,  and  calculating  control  gains  (see  [1]).  This  is  accomplished  by  approximating 
the  beam  displacement  and  acoustic  potential  by  spline  and  spectral  expansions,  respectively. 

3.1  System  Approximation 

As  detailed  in  [1]  where  the  corresponding  linear  problem  Weis  considered,  crvbic  splines  are 
suitable  for  discretizing  the  beam  displacement  since  they  satisfy  the  smoothness  requirement 
as  well  as  being  easily  implemented  when  adapting  to  the  fixed-end  boundary  conditions  and 


patch  discretizations.  Specifically,  the  approximate  beam  displacement  is  taken  to  be  the 
linear  combination 

i=l 

where  is  the  cubic  spline  which  has  been  modified  to  satisfy  the  boundary  conditions. 
The  acoustic  potential  is  approximated  by  the  Galerkin  expansion 


j=0  t=o 
i+J#0 

where  Pi{x)  and  P/(j/)  denote  the  standard  Legendre  polynomials  that  have  been  scaled 
by  transformation  to  the  intervals  [0,a]  and  [0,^],  respectively.  The  condition  i  +  j  /  0 
eliminates  the  constant  function  thus  guaranteeing  that  the  set  of  functions  is  suitable  as  a 
basis  for  the  quotient  space.  We  take  {Bk{x,^)}f-i  =  where  m  = 

(mi  -f  l)(my  +  1)  -  1. 

The  m  and  n  —  1  dimensional  approximating  cavity  and  beam  subspaces  are  taken  to  be 
=  span{5f‘})^j  and  =  span{B"}”rj\  respectively,  where  By  and  Bf  are  the 
cavity  and  beam  bases  described  above.  Defining  N  =  m  +  n  —  I,  the  approximating  state 
space  is  then  taken  to  be  =  Hy  x  and  the  product  space  for  the  first  order  system  is 
=  X  H^. 

By  restricting  the  infinite  dimensional  system  (2.11)  to  xH^ ^  one  obtains  the  nonlinear 

finite  dimensional  system 


■  My 

0  1 

■  d^{t)  ■ 

0 

My 

. 

■  My 

0  ] 

■  t?^(o)  ■ 

0 

My  \ 

_ 

0  uy  1  r 

-Ay  -Ay 


u{t)  + 


0 

F^{t) 


My  =  diag[M^[,  My]  , 


=  diagM,  MS  , 


,  F''(i)  =  [0,Fi'(()]'  . 

The  vector  (f>2  {t)i  ^^(0?  ^^(0?  ‘  ‘  ‘  i  ^n^-i  i^)V  contains  the  Nxl  = 

(m  +  n  —  1)  X  1  approximate  state  coefficients  while  u(<)  =  (ui(t),  •  •  •  Us(<))  contains  the  .s 
voltage  values.  The  matrices  M^  and  are  the  mass  and  stiffness  matrices  which  arise  when 
solving  the  uncoupled  wave  equation  with  Neumann  boundary  conditions  while  M^,  and 
A22  the  mciss,  stiffness  and  damping  matrices  which  arise  when  solving  the  damped  beam 
equation  with  fixed  boundary  conditions.  The  matrices  M^  and  M^  result  from  the  choice 
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of  V  inner  product  (see  [1]).  The  contributions  from  the  coupling  terms  are  contained  in  the 
matrix  and  operator  while  the  control,  forcing  and  initial  terms  are  contained 

in  B2  1  F-iit)  and  respectively.  A  more  detail  description  of  the  various  component 
matrices  can  be  found  in  [1]. 

The  equivalent  finite  dimensional  Cauchy  equation  is  then 

yN{t)  =  .4^  [yN[t))  +  u{t)  +  {t)  ^ 

with  the  various  operator  definitions  following  from  those  given  above.  As  with  the  infinite 
dimensional  problem,  it  will  also  be  useful  to  consider  the  corresponding  problem 

=  AVit)  +  B^u{t)  +  F^{t) 

y^(0)  =  Jo" 

when  considering  issues  such  as  parameter  convergence.  The  linear  operator  A^  is  obtained 
by  considering  the  linearization  of  cis  discussed  in  [1]. 

3.2  Algorithm  for  Constructing  the  Nonlinear  Component 

The  determination  of  state  trajectories  involves  the  repeated  construction  of  the  nonlinear 
operator  Although  most  of  the  components  of  this  operator  need  to  be  constructed 

only  once,  the  (n  —  1)  x  m  matrix  niust  be  determined  at  each  step  in  the  solution 

of  the  ordinary  differential  equation  (ODE)  (3.1)  and  hence  its  formation  must  be  made  as 
efficient  as  possible. 

We  first  point  cut  that  for  j  =  1,  •  •  • ,  mj,,  rriy,  this  (n  —  1)  x  m  matrix  heis  the  entries 

[<  K(*))]„..  =  /„“  B;{x)Pnx)Pl  (e  w^{t)B'}{x)j  dx 
where  the  integrals  are  evaluated  via  a  Gaussian  quadrature  rule  of  the  form 

n, 

f{x)dx  ft!  Y,  Ckf{xk)  . 

k=l 

Here  n,  is  the  number  of  quadrature  points  and  Ck  is  the  quadrature  weight.  These  quadratures 
and  hence  the  formation  of  the  matrix  can  be  efficiently  accomplished  in  the  manner  outlined 
in  the  following  algorithm. 

Numerical  Algorithm  for  Creating  the  Nonlinear  Component  (w^  (')): 
i.  Create  the  matrices  ,  A^2a'w  ^Z2bj\  which  have  the  components 

(hence dim  =  n,x(n-l)  ,  dim  A^2aw  =  («-l)xng  and  dim  A^2hj\  =  n,x(ma:+l)). 
These  matrices  need  to  be  formed  only  once  and  can  be  created  before  solving  the  ODE 


2.  With  the  notation  o  denoting  Hadamard  or  componentwise  multiplication,  perform  the 
following  operations  in  the  ODE  solver: 

X  =  X  =  7lg  X  1 

for  j  =  1  :  niy  +  1 

b  =  Pjix)  ,  dim  b  =  iiq  X  \ 

^32bj-2  =b*  07ies{l,7na;  +1)  ,  dim  X  +  1) 

^32bj  =  ^326jl  °  ^^^26^2  (^^(O)  >  dim  A^^bj  X  {771^  +  1) 

mdl  =  (i  -  1)  •  (mi  +  1)  +  1 
i7id2  =  j  •  (rui  +  1) 

{•,indl  :  md'l)  =  A^2bj  »  dim  A^^b  («'^(0)  =  x  7n 

end 

(«'^(*))  =  Pf^32aw  •  >i^25  >  dim  >1^2  =  (n  -  1)  X  m  . 

By  creating  the  matrices  ,  A^2aw  A^2bj\  offline,  the  time  needed  to  solve  the  ODE 
system  is  reduced  thus  improving  the  efficiency  of  the  scheme  for  parameter  estimation  and 
control  applications. 

4  Parameter  Estimation  Problem 

The  goal  of  the  parameter  estimation  problem  is  to  determine  the  “true”  material  parameters 
pfc,  jB/,  cdI  and  /Cf ,  /Cf  ,  i  =  1,  •  •  • ,  s,  given  data  meaisurements  z  from  some  observable  sub¬ 
space  Z  of  the  state  space.  To  pose  this  mathematically,  we  let  q  =  (pt,  EI,CDl,fCf,  •  •  •  ,X!if , 
/Cf , •  •  • , JCg)  and  assume  that  q  €  Q  where  Q  denotes  an  appropriately  chosen  admissible 
parameter  space.  The  parameter  estimation  problem  is  to  then  seek  q  E  Q  which  minimizes 

J(,)  =  |c,[c,{^(i,,,)}-l.]|"  (4.3) 

given  pointwise  temporal  meeisurements  z,  =  z{ti)  at  given  points  on  the  beam  and  in  the 
cavity,  Note  that  this  minimization  is  performed  subject  to  Z  =  {<j>,x,^,7j})  satisfying  the 
coupled  system  equations  (2.11)  or  (2.15).  Depending  on  the  experimental  apparatus,  the 
data  observations  may  consist  of  position,  velocity,  acceleration,  or  accumulated  strain  mea¬ 
surements  at  points  on  the  beam  as  well  as  pressure  me£isurements  inside  the  cavity.  The  form 
of.  the  operator  C2  depends  upon  whether  one  is  performing  the  estimation  procedures  in  the 
time  domain  or  in  the  frequency  domain.  In  the  time  domain,  C2  is  the  identity  whereas  it  is 
the  Rsurier  transform  for  estimation  in  the  frequency  domain. 

For  time  domain  estimation  with  data  consisting  of  position,  velocity,  or  acceleration 
nieasurements  at  points  x  on  the  beam,  the  fit  criterion  functional  to  be  minimized  is 

2 

•^(</)  =  H 
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with  1/  =  0, 1  or  2.  On  the  other  hand,  if  a  patch  with  endpoints  at  x\  and  X2  is  used  to  collect 
accumulated  strain  measurements,  an  appropriate  functional  is 


J(i)  =  E 


{dx 


dw 

dx 


(4.5) 


(see  (2.9)).  Here  the  data  consists  of  the  voltage  measured  across  the  patches.  The  patch  data 
can  also  be  combined  with  pressure  measurements  p,  at  points  {x,  y)  in  the  cavity  to  provide 
a  fit  criterion  functional 


J{<!)  =  E 


4-  \<f>t{U,x,y)-pi 


(4.6) 


which  in  some  cases  has  more  sensitivity  than  that  in  (4.5)  which  considers  only  strain  mea¬ 
surements. 

The  above  fit  criteria  can  also  be  used  with  data  that  has  been  transformed  to  the  frequency 
domain  (in  which  case  C2  is  the  Fourier  transform),  and  this  is  indeed  a  common  procedure 
for  starting  the  optimization  process  with  data  in  which  several  frequencies  are  excited  (see 
the  comments  in  the  examples  as  well  as  [13,  14]).  In  this  case,  optimization  is  qualitatively 
performed  in  the  frequency  domain  until  frequencies  match  sufficiently  so  that  the  optimization 
routines  will  converge  with  the  time  domain  data. 

To  facilitate  the  estimation  of  the  material  parameters  p,  El  and  cp/,  we  now  make  some 
assumptions  regarding  their  spatial  behavior.  Because  the  beam  and  patches  are  considered 
to  be  homogeneous  as  well  as  uniform  in  width  and  thickness,  it  is  reasonable  to  assume  that 
the  density,  stiffness  and  damping  parameters  of  the  combined  beam/piezoceramic  patches  are 
piecewise  constant  in  nature  (see  for  example,  [15]).  A  suitable  partition  is  then  taken  to  be 
{xk}  =  {0,a}  where  the  2s  points  {a,j}  are  the  endpoints  of  the  s  piezoceramic 

patches.  Finally,  we  assume  that  these  parameters  have  the  form 


p{x) 

2s+l 

=  E  CkBk{x) 

k=l 

II 

II 

0 

’  •  =  C2s+i 

EI{x) 

2s+l 

=  E)  ^kBk{x) 

fc=l 

II 

II 

■  •  =  C2S+1 

cdI{x) 

2s+l 

=  E  ^kBk{x) 

fc=l 

,  Cl  =  C3  =  • 

■  '  =  C2s+1 

where  the  piecewise  constant  basis  functions  are  defined  by  Bk{x)  =  H{x  —  Xk-i)  —  H{x  —  Xk). 
The  coefficient  constraints  cj  =  C3  =  •  •  •  =  C2s+i,  and  so  on,  result  from  the  uniformity  of  the 
beam  in  areas  not  covered  by  patches.  Finally,  we  recall  from  the  definitions  (2.8)  and  (2.9) 
that  the  patch  parameters  /Cf , •  •  •  ,/Cf ,  and  /Cf, •  •  •  ,/Cf  are  simply  constants  which  depend 
on  piezoelectric  properties,  the  geometry  and  size  of  the  patch,  and  patch  and  bonding  layer 
properties. 

Although  Q  is  finite  dimensional  with  the  above  assumptions  on  the  parameters,  the  min¬ 
imization  of  the  fit  criteria  in  (4.4),  (4.5)  and  (4.6)  involves  an  infinite  dimensional  state 
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and  hence  is  not  immediately  tractable.  With  the  Galerkin  schemes  of  the  the  last  section, 
however,  corresponding  minimization  problems  involving  the  state  approximations  can  be  de¬ 
veloped  and  used  when  estimating  physical  parameters  with  these  fit-to-data  techniques.  With 
<f>^  and  C  V  defined  in  the  last  section,  finite  dimensional  functionals  corresponding 
to  those  in  (4.4),  (4.5)  and  (4.6)  are 


=  E 

J%)  =  E 


XV 


df' 


(t,',  Xj  (j)  Zx 


,  7/  =  0, 1,2  , 


.c  ( dtv^  dxv^ .  \ 

(^15^2)?)  [til  </)  1 


\  dx 


and 

J 


div^ 


dxv^ 


{U, X2;  q)  -  <7)  1  - 


+  \(l>^{U,Xxy)-pi\^\  , 


(4.8) 

(4.9) 


(4.10) 


respectively.  The  approximate  beam  displacement  10^  and  acoustic  pressure  pf<j>^  at  the 
various  points  are  found  by  solving  either  the  nonlinear  finite  dimensional  system  (3.1)  or  the 
system  (3.2)  if  one  is  considering  the  linearized  problem. 

The  following  theorem  taken  from  [16]  specifies  conditions  under  which  convergence  and 
continuous  dependence  (on  data)  of  the  solutions  to  the  linearized  finite  dimensional  parameter 
estimation  problems  involving  the  functionals  (4.8),  (4.9)  and  (4.10)  can  be  expected. 


Theorem  2.  (Linearized  System)  Let  (5  be  a  a  compact  subset  of  a  metric  space  Q  with 
metric  d  and  assume  that  C  V  approximates  V  in  the  sense  that  for  each  $  G  K, 
there  exists  G  such  that 


N 


<  £{N)  0  as  N  ^  00  . 


(4.11) 


Furthermore,  assume  that  (Xi{q)  and  <T2(^)  defined  in  (2.12)  are  K-elliptic,  continuous, 
and  satisfy  the  continuity  with  respect  to  parameter  condition 


lo-,(^)(<>,'I')  -  (ri(7)($,^)|  <  7.%,^)|$|v|^|v  ,  for  G  V 
for  i  =  1,2  and  q,q  E  Q.  Finally,  assume  that 

q  i-*  {Bu  -f  F){t;q)  is  continuous  from  Q  to  L^((0,!r),  K*)  . 
For  arbitrary  such  that  q’^  —*  q  in  Q,  one  then  has  the  convergence 

z^{t-,q^)  — >  z{t\q)  in  V  norm 

»  ZiiUq)  in  y  norm 


(4.12) 


(4.13) 


(4.14) 


for  <  >  0.  Here  z  and  Zt  are  solutions  to  the  linearized  system  (2.19)  and  z^  and  z^ 
solve  the  corresponding  linear  finite  dimensional  system  in  . 
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In  our  problems  the  admissible  parameter,  space  Q  is  taken  to  be  a  compact  subset  of 
the  metric  space  Q  =  [Loo(0,a)]^  x  [IR*]^  with  elements  satisfying  the  coefficient  constraints 
specified  in  (4.7)  as  well  as  the  physical  constraints  p  >  0  ,  El  >  0  and  cdI  >  0  on  (0,  o).  In 
considering  the  remaining  hypotheses,  we  note  that  the  sesquilinear  forms  satisfy  the  ellipticity, 
continuity  and  parameter  continuity  conditions  as  long  as  damping  is  assumed  in  some  region 
n  C  n  (see  (2.13)  and  (2.14)).  Moreover,  the  input  term  (F  +  Bu){q)  satisfies  the  condition 
(4.13).  Finally,  the  convergence  condition  specified  in  (4.11)  is  satisfied  as  a  consequence  of 
the  approximating  properties  of  the  cubic  splines  and  Legendre  polynomials  in  a  Galerkin 
setting  (see  [6]  for  further  details). 

Hence,  for  the  linearized  problem  with  acoustic  damping  in  fl,  a  subsequence  of  solutions 

to  the  problems  involving  the  minimization  of  the  functionals  (4.8)  with  //'  =  0, 1 ,  (4.9)  or 
(4.10)  subject  to  (3.2)  will  converge  to  a  solution  q  of  the  original  problem  of  minimizing  the 
functionals  (4.4),  (4.5)  or  (4.6)  subject  to  (2.19).  The  convergence  in  the  case  involving  the 
minimization  of  (4.8)  with  acceleration  data  does  not  follow  directly  from  this  theorem  but 
can  be  obtained  using  results  from  the  theory  of  analytic  semigroups  in  a  manner  analogous 
to  that  used  in  [7]. 

With  boundedness  and  Lipschitz  continuity  assumptions  on  the  nonlinear  coupling  term 
G{Z),  similar  results  can  be  obtained  for  the  nonlinear  problem  as  summarized  in  the  following 
remark. 

Lemma  2.  (Nonlinear  Damped  System)  Consider  the  system  with  the  nonlinear  input 
term  Bu{t)  +  F{t)  +  G{Z{t)).  If,  in  addition  to  assuming  a  continuity  condition  of  the 
form  (4.13)  (with  Bu  +  F  replaced  by  BU+  F  +  G{Z)),  we  also  assume  that  G{Z{t))  is 
continuous  in  t,  uniformly  Lipschitz  continuous  in  Z,  and  displays  at  most  affine  growth 
at  oo,  then  convergence  results  analogous  to  those  summarized  in  Theorem  2  can  be 
obtained  for  the  nonlinear  system. 

We  point  out  that  these  continuity  assumptions  on  G  were  also  made  when  discussing  the  well- 
posedness  of  the  nonlinear  system.  Details  concerning  these  conditions  as  well  as  arguments 
leading  to  the  proof  for  the  nonlinear  case  can  be  found  in  [5]. 

As  indicated  previously,  the  assumption  of  medium  damping  inside  the  cavity  is  often 
inappropriate  in  applications  of  interest.  While  we  have  not  extended  Theorem  2  to  include 
the  case  in  which  acoustic  damping  is  omitted,  extensive  numerical  tests  have  indicated  that 
parameter  convergence  and  continuous  dependence  of  the  parameters  on  data  is  being  obtained 
in  the  same  manner  exhibited  by  the  system  having  both  acoustic  and  structural  damping. 
This  is  demonstrated  by  the  results  in  the  following  examples  for  the  nonlinear  2-D  structural 
acoustic  system  in  which  cavity  damping  is  omitted. 
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5  Numerical  Examples 

To  test  the  parameter  estimation  methodology,  the  general  problem 
(f)ti  =  c^A(f)  (a:,j/)  6  n  >  0  , 

'71  =  0  (x,y)  E  r  ,t  >  0  . 


V<f>{t,x,w{t^x))  •  n  =  pfWt{t,x)  0  <  X  <  .0  >  0  , 


pwit  + 


f„rd^w 


EI^  +  cdI 


dx^\  dx^ 


dx^dt 


-<f)t{t,x,w{t,x))  +  f{t,x)  +  ■^j2K:fui{t)xpei{x) 


w{t,0)  = -^{t,0)  =  w{t,  .0)  = -^{t,  .6)  =  0  t>0  , 


<t>{0,x,y)  =  ij>t{0,x,y)  =  txj(0,x)  =  iU4(0,a;)  =  0 


0  <  X  <  .6  , 

^>0, 


was  considered.  The  cavity  was  assumed  to  have  x  and  y  dimensions  a  =  .6  m  and  £  =  \  m 
with  a  beam  at  one  end  having  length  .6  m,  width  .1  m  and  thickness  .005  m  (see  Figure  3). 
The  density  and  Young’s  modulus  for  regions  of  the  beam  devoid  of  patches  were  taken  to  be 
Pb  —  2700  kglm^  and  £;  =  7.1  x  10^“  iV/m^  which  yields  p  =  1.35  kgim  and  El  =  73.96  Nni^ 
for  the  linear  mass  density  and  stiffness  parameter  (see  Table  1  for  a  compilation  of  the 
structural  parameters  for  the  system).  The  Kelvin- Voigt  damping  parameter  was  chosen  to 
be  cpI  =  .001  kgm^fsec.  Finally,  the  values  pj  =  1.21  kglm^  and  c  =  343  mfsec  were  used 
for  the  atmospheric  density  and  speed  of  sound. 

In  the  examples,  we  consider  a  system  in  which  the  bounding  end  beam  hais  bonded 
to  it  a  centered  piezoceramic  patch  covering  1/3  of  its  length  as  shown  in  Figure  3.  The 
patch  is  cissumed  to  have  thickness  T  =  .000508  m  and  width  6  =  .1  m  (we  point  out 
that  the  chosen  thickness  value  corresponds  to  20  mil  which  is  a  commercially  available 
thickness  for  piezoceramic  patches).  The  Young’s  modulus  and  density  were  taken  to  be 
Epe  =  6.3  X  10^°  Nj 771?  and  ppe  =  7650  kgfm^  which  are  reasonable  for  a  patch  made  from 
G-1195  piezoceramic  material. 

From  (2.7)  and  (2.8),  we  see  that  the  density  and  stiffness  coefficient  in  the  region  of  the 
combined  beam  and  patch  (Region  2)  will  be  greater  than  that  of  the  beam  (Region  1)  (see 
Figure  3).  We  also  assume  that  the  damping  coefficient  will  be  slightly  larger  in  Region  2 
than  Region  1. 

For  testing  purposes,  the  structural  parameter  values  in  regions  covered  by  the  patches  were 
chosen  as  specified  in  Table  1.  As  seen  there,  the  constant  /C®  =  Epebd3\{h  +  2Tbe  +  T),  which 
arises  when  modeling  the  actuation  due  to  the  patch,  Wcis  taken  to  be  JC^  =  .0067  NmfV 
(this  latter  value  was  obtained  by  assuming  a  bonding  layer  of  thickness  Tbe  =  .0001  m  and 
taking  d^i  =  1.9  x  10“^°  ni/V  which  is  the  value  specified  for  G-1195  piezoceramic  material). 
The  constant  IC^  which  arises  when  using  the  patches  as  sensors  weis  taken  to  have  the  value 
=  170  V. 


22 


Fig.  3.  Acoustic  chamber  with  one  centered  1/3  length  patch. 

The  following  examples  demonstrate  the  numerical  estimation  of  the  material  parameters 
p,  iS/jCc/j/C®  and  JC^  using  various  techniques  for  exciting  the  system  and  observing  the 
response.  In  the  first  example,  the  naturaLirequencies  for  the  fully  coupled  system  were  de¬ 
termined  by  simulating  an  impulse  hammer  impact  to  the  center  of  the  beam.  The  knowledge 
of  these  natural  frequencies  was  then  used  when  choosing  the  frequencies  of  the  exciting  forces 
in  the  other  examples.  In  the  second  example,  a  periodic  uniform  forcing  function  (model¬ 
ing  a  uniform  exterior  acoustic  pressure  field)  was  applied  to  the  beam  for  a  short  interval 
of  time  and  then  set  to  zero.  This  forcing  function  was  chosen  so  that  three  system  modes 
were  initially  excited  and  then  allowed  to  begin  decaying  due  to  the  damping  the  beam.  The 
acceleration  of  the  center  of  the  beam  was  used  as  data  for  estimating  the  material  parameters 
pj  El  and  cdI  (since  no  voltage  was  applied  in  this  example  and  the  patch  was  not  used  for 
sensing,  and  K,^  were  not  estimated).  In  the  third  example,  the  system  was  excited  by 
the  application  of  a  periodic  voltage  into  the  patch.  Again,  the  system  was  excited  for  a  short 
time  interval  and  then  allowed  to  freely  decay  in  energy.  Acceleration  data  at  the  center  of 
the  beam  was  used  to  estimate  the  four  parameters  p,  £/,  c©/  and  /C®.  The  patch  was  used 
both  for  actuating  and  sensing  in  the  fourth  and  fifth  examples.  In  Example  4,  a  periodic 
voltage  was  applied  for  an  initial  time  interval  after  which  the  system  energy  was  allowed  to 
decay.  During  the  decay  interval,  two  sets  of  data  were  calculated  and  a  comparison  was  made 
between  the  results  obtained  when  each  was  used  for  recovering  the  parameters  p,  El,  cq!,  K.^ 
and  K^.  The  first  data  set  consisted  solely  of  the  voltage  produced  by  the  patches  during  vi¬ 
bration  (and  hence  contained  strain  measurements)  while  the  second  contained  a  combination 
of  voltage  measurements  from  the  patch  and  acoustic  pressure  values  from  inside  the  cavity. 
Finally,  a  simulated  voltage  spike  to  the  patches  was  used  to  excite  the  system  in  Example  5 
(with  an  effect  similar  to  that  observed  when  an  impulse  hammer  is  used  to  excite  the  system) 
with  the  patches  again  being  used  as  a  sensor  throughout  the  remainder  of  the  time  interval. 
Thus  in  the  last  two  examples,  the  “smart  material”  aspects  of  the  structure  were  utilized  in 
determining  its  physical  parameters. 
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True  Values 

Initial  Values 
(Set  1) 

Initial  Values 
(Set  2) 

p 

Region  1 

1.35 

1.4 

1.4 

(%/m) 

Region  2 

2.115 

2.0 

2.0 

El 

Region  1 

73.96 

74.0 

75.0 

{Nm?) 

Region  2 

125.4 

125.0 

127.0 

col 

Region  1 

.001 

.0008 

.0001 

{Nm'^sec) 

Region  2 

.00125 

.0008 

.0001 

JC^  {N 

mfV) 

.0067 

.007 

.01 

JC^{V) 

170.0 

172.0 

175.0 

Table  1.  (a)  True  values  of  the  material  parameters;  (b)  Set  1  of  initial  guesses  for  the  material 
parameters;  (c)  Set  2  of  initial  guesses  for  the  material  parameters. 


Examjple  1:  System  Dynamics 


in  order  to  determine  the  system  dynamics  with  the  parameter  values  in  Table  1,  the 
forcing  function  /  was  chosen  to  simulate  an  impuls  *  at  the  center  of  the  beam;  that  is, 


f{t,x)  =  S{x-.Z)6{t). 

This  models  the  force  that  would  be  delivered  by  a  centered  impulse  hammer  hit.  After  the 
initial  impulse,  the  system  was  allowed  to  run  unforced  through  time  T  =  8/60. 

The  beam  acceleration  obtained  with  rux  =  my  =  12  and  n  =  16  basis  functions  at  the 
point  X  =  .3  is  plotted  in  Figure  4  with  a  corresponding  frequency  plot  in  Figure  5.  The  first 
four  system  responses  occur  at  62.9, 179.7,342.1  and  397.9  hertz. 

For  comparison,  we  note  that  the  analytic  natural  frequencies  of  the  first  two  symmetric 
modes  of  an  isolated,  homogeneous  (no  patches),  undamped  beam  having  the  same  dimensions 
as  those  in  this  system  are  73.2  and  395.6  hertz  while  those  of  the  uncoupled  acoustic  cavity 
are  171.5  and  343  hertz  (the  analytic  undamped  beam  and  cavity  frequencies  are  given  by 


/.= 


A? 

2^0^  V  P 


Ai  =  4.7300,  A2  =  7.8532,  •  •  • 


(5.1) 


and 


respectively).  The  differences  between  the  observed  and  analytic  values  are  due  to  the  presence 
of  the  patches  on  the  beam,  the  internal  damping  in  the  beam,  and  the  coupling  between  the 
beain  and  acoustic  field  dynamics. 
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The  frequencies  observed  here  can  also  be  compared  with  the  system  values  65.9, 181.3, 
343.9  and  387.8  hertz  which  were  obtained  when  no  patches  were  bonded  to  the  plate  and 
linearized  coupling  conditions  were  assumed  (see  [8]).  It  can  be  seen  that  the  increased  stilfness 
due  to  the  presence  of  the  patches  manifests  itself  at  the  higher  frequencies  (397.9  hertz  versus 
387.8  hertz  for  the  uniform  beam)  whereas  the  increased  density  is  more  of  a  factor  at  the 
lower  frequencies  (62.9  hertz  versus  65.9  hertz  for  the  uniform  beam).  This  can  be  explained 
by  comparing  the  bending  shapes  of  the  first  and  third  beam  modes  and  noting  that  the 
increased  stiffness  is  more  of  an  influence  on  the  third  beam  mode  while  the  increased  density 
more  directly  influences  the  first  mode. 

Beam  Response  to  an  Initial  Impulse 


lime  (sec) 


Fig.  4.  The  beam  response  to  a  centered  impact. 


Beam  Response  to  an  Initial  Impulse 


Frequency  (hertt) 


Fig.  5.  The  beam  frequency  response  to  a  centered  impact. 
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Example  2:  Periodic  Acoustic  Excitation,  Acceleration  Data 


In  this  example,  a  uniform  (in  space)  periodic  force  modeling  an  acoustic  plane  wave  was 
used  to  excite  the  system  and  acceleration  data  was  used  to  estimate  the  parameters  p,  El 
and  C£)/  (/C®  and  JC^  were  not  estimated  since  we  were  not  applying  voltage  to  the  patch  or 
using  the  patch  for  sensing).  Specifically,  the  forcing  function  was  taken  to  be 


sin(1207rf)  +  sin(3607r<)  +  sin(8007rt)  ,  0  <  t  <  1/60 

0  ,  1/60  <8/60 


which  initially  excites  the  first,second  and  fourth  system  modes  and  then  allows  the  oscillations 
to  begin  dying  away  due  to  the  damping  in  the  beam  (see  Figures  6  and  7). 

The  parameter  estimation  was  performed  with  data  which  was  generated  by  calculating  the 
acceleration  of  the  central  point  of  the  beam  at  498  uniformly  distributed  points  throughout 
the  time  interval  [0,8/60]  (hence  i/  =  2  in  (4.8)).  The  acceleration  was  determined  by  using 
a  second-order  central  difference  on  the  displacements  which  were  obtained  by  solving  the 
nonlinear  finite  dimensional  system  (3.1)  and  evaluating  at  the  points  (.3,4),  4  = 
k-  eo(5oo)’  ^  =  2,  •  •  •  ,499.  Due  to  the  relatively  small  number  of  frequencies  being  matched,  all 
identification  procedures  were  performed  in  the  time  domain  which  implies  that  C2  in  (4.3) 
was  taken  to  be  the  identity. 

To  test  the  algorithm  and  software,  noisefree  data  was  first  generated  using  120  acoustic 
and  15  beam  basis  functions  =  10,  n  =  16).  Using  the  initial  parameter  choices 

in  column  5  of  Table  1,  the  parameters  were  estimated  using  the  same  basis  choice  with  the 
optimization  being  performed  via  a  Levenberg-Marquardt  routine.  As  dem  .'‘strated  by  the 
recovered  values  in  the  third  column  of  Table  2,  very  accurate  estimates  ot  the  parameters 
can  be  obtained  when  using  the  same  number  of  basis  functions  when  generating  the  data 
and  estimating  the  physical  parameters. 

To  provide  more  realistic  simulations,  data  was  then  generated  using  168  acoustic  and 
19  beam  basis  functions  {rrix  =  mj,  =  12,  n  =  20)  with  120  acoustic  and  15  beam  basis 
functions  again  being  used  for  estimating  the  parameters  (the  use  of  a  larger  number  of  basis 
functions  when  generating  the  data  has  the  effect  of  adding  numerical  noise  to  the  values  being 
approximated  in  the  optimization  routine).  The  parameter  estimates  obtained  with  noisefree 
data  and  data  to  which  10%  relative  noise  was  added  are  reported  in  Table  2.  While  these 
results  were  obtained  with  the  initial  values  in  column  5  of  Table  1,  they  are  representative 
of  those  obtained  with  a  variety  of  initial  guesses  (for  this  method  of  system  excitation  and 
observation,  the  optimization  routine  converged  for  a  relatively  large  range  of  initial  values). 
Time  and  frequency  domain  plots  of  the  data  containing  10%  noise  and  the  acceleration 
obtained  with  the  estimated  parameter  values  of  Table  2,  column  5,  are  given  in  Figure  6 
and  7.  We  note  that  in  these  figures  it  is  essentially  impossible  to  distinguish  between  the 
time  data  and  the  model  response  with  estimated  parameters  with  the  graphical  resolution 
used.  These  plots  reinforce  the  observation  that  reaisonable  estimates  of  the  parameters  can  be 
obtained  with  acoustic  excitation  and  acceleration  data  with  or  without  noise  in  the  simulated 
data. 
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No  Noise 
(Data  Set  1) 

No  Noise 
(Data  Set  2) 

10%  Noise 
(Data  Set  2) 

p 

Region  1 

1.3500 

1.3499 

1.3469 

{kglm) 

Region  2 

2.1150 

2.1149 

2.1181 

El 

Region  1 

73.960 

73.836 

72.130 

[Nm^) 

Region  2 

125.400 

125.584 

130.716 

CD  I 

Region  1 

.001000 

.001005 

.001023 

{Nm^sec) 

Region  2 

.001250 

.001232 

.001091 

Table  2.  Estimated  values  of  p,  El  and  cdI;  (a)  Data  generated  with  nix  =  7ny  =  10, 7i  =  16  and 
no  noise  added  to  data;  (b)  Data  generated  with  iiix  =  niy  —  \2,n  =  20  and  no  noise  added  to  data; 
(c)  Data  generated  with  nix  =  =  12,  ti  =  20  and  10%  noise  added  to  the  data. 


System  Response  to  an  Initial  Periodic  Force 


Fig.  6.  The  system  response  to  an  initial  periodic  force.  The  solid  line  denotes  the  data  containii  .g 
10%  noise  while  the  dashed  line  is  the  model  response  obtained  with  the  reeovered  parameter  values 
in  column  5  of  Table  2. 
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System  Response  to  an  Initial  Periodic  Force 


Fig.  7.  The  system  frequency  response  to  an  initial  periodic  force.  The  solid  line  (with  x ’s)  denotes 
the  data  containing  10%  noise  while  the  dashed  line  (with  o’s)  is  the  model  response  obtained  with 
the  recovered  parameter  values  in  column  5  of  Table  2. 


Example  3:  Periodic  Patch  Excitation,  Acceleration  Data 


A  second  means  of  exciting  the  system  is  through  the  application  of  an  out-of-phase  voltage 
to  the  piezoceramic  patches.  The  driving  voltage  here  was  taken  to  be 


u{t)  = 


sin(1207ri)  sin(3607rt)  -|-  (.l)sin(8007rt) 

0  ,  1/60  <  t  <  8/60 


0  <  t  <  1/60 


which  initially  excites  the  first  four  system  modes  and  then  allows  the  oscillations  to  begin 
dying  away  due  to  the  damping  in  the  beam  (see  Figures  8  and  9).  Because  the  patch  excites 
the  higher  frequencies  more  efficiently  than  did  the  simulated  sound  field  of  the  last  example, 
the  magnitude  of  the  400  hertz  contribution  was  reduced  to  better  balance  the  system  response. 

Acceleration  data  obtained  in  the  same  manner  described  in  the  last  example  was  used  to 
estimate  the  four  parameters  p,  El,  cdI  and  /C®  and  the  estimates  obtained  with  no  noise  and 
10%  noise  added  to  the  data  are  summarized  in  Table  3  (see  Table  1  for  the  true  parameter 
values).  Figures  8  and  9  contain  the  time  and  frequency  domain  results  for  the  case  in 
which  the  data  contains  10%  noise.  In  both  cases,  the  choices  rtix  =  rriy  =  12,n  =  20 
and  rtix  =  rtiy  —  10,  n  =  16  were  used  when  generating  the  data  and  estimating  parameters, 
respectively.  Also,  while  similar  results  were  obtained  with  a  relatively  large  variety  of  starting 
parameter  values,  the  parameter  values  reported  here  were  obtained  with  the  initial  choices 
in  column  3  of  Table  1. 


No  Noise 

10%  Noise 

p 

Region  1 

1.3502 

1.3425 

{kgfm) 

Region  2 

2.1211 

2.0973 

El 

Region  1 

73.139 

73.449 

(Nm^) 

Region  2 

128.179 

124.741 

CD  I 

Region  1 

.000997 

.000849 

(Nm^sec) 

Region  2 

.001256 

.001787 

K.^{N 

m/V) 

.006771 

.006654 

Table  3.  Estimated  values  of  p,EI,cdI  and  fC^;  (a)  No  noise  added  to  the  data;  (b)  10%  noise 
added  to  the  data. 

In  comparing  the  results  obtained  with  noisefree  data  with  those  obtained  under  the  same 
conditions  (same  number  of  basis  functions  and  noisefree  data)  in  the  Icist  example,  we  note 
an  approximately  2%  change  in  the  estimated  value  of  the  stiffness  parameter  in  the  region 
covered  by  the  patches.  This  difference  appears  to  be  due  to  the  fact  that  the  “true”  data 
is  calculated  with  a  larger  number  of  basis  functions  and  hence  with  more  accuracy  than  the 
solutions  being  obtained  during  the  parameter  estimation.  The  resulting  “numerical  noise” 
manifests  itself  more  strongly  in  this  case  since  the  patches  are  more  effective  than  the  acoustic 
field  at  exciting  high  frequency  oscillations  which  require  greater  accuracy  to  resolve  (this 
tendency  is  also  noted  in  the.  remaining  examples  where  patch  activation  is  used  to  excite  the 
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system).  We  point  out  that  when  the  same  number  of  basis  functions  are  used  for  generating 
data  and  estimating  parameters,  the  results  are  essentially  identical  to  those  in  column  2  of 
Table  2  with  highly  accurate  estimates  of  the  physical  parameters.  Finally  we  note  that  the 
“numerical  noise”  due  to  the  differing  discretizations  combines  with  random  noise  added  to  the 
data  to  simulate  the  noise  which  is  present  when  real  data  is  used  to  estimate  the  parameters 
in  a  physical  experiment. 

\ 

System  Response  to  an  Initial  Periodic  Voltage 


002  0  04  0  06  0.08  0.1  012 

Time  (sec) 


Fig.  8.  The  system  response  to  an  initial  periodic  voltage.  The  solid  line  denotes  the  data  containing 
10%  noise  while  the  dashed  line  is  the  model  response  obtained  with  the  recovered  parameter  values 
in  Table  3. 


System  Response  to  an  Initial  Periodic  Voltage 


Frequency  (hertz) 


FiG;  9.  The  system  frequeticy  response  to  an  initial  periodic  voltage.  The  solid  line  (with  x's) 
denotes  the  data  containing  10%  noise  while  the  dashed  line  (with  o’s)  is  the  model  response  obtained 
withlthe' recovered  parameter  values  in  Table  3. 


Example  4:  Periodic  Patch  Excitation,  Patch  and  Pressure  Data 


In  the  previous  two  examples,  acceleration  data  modeling  that  which  would  be  obtained 
from  a  centered  accelerometer  was  used  in  the  estimation  of  the  physical  parameters.  In  this 
as  well  cis  the  next  example,  both  the  excitation  of  the  system  and  the  sensing  of  the  beam 
dynamics  are  performed  with  the  patches,  thus  utilizing  their  “smart  material”  capabilities. 
The  beam  data  in  this  case  consists  of  the  voltage  that  is  produced  when  the  patches  strain 
during  vibration  (hence  the  criterion  functional  (4.9))  as  well  as  data  consisting  of  voltage 
measurements  in  combination  with  pressure  measurements  from  within  the  cavity  (with  the 
criterion  functional  (4.10)). 

Here  a  periodic  driving  voltage  was  applied  to  the  patches  for  the  first  1/5  of  the  time 
interval  [0, 10/60]  of  interest  after  which  the  system  wa.  allowed  to  begin  decaying  in  energy. 
Specifically,  the  input  voltage  was  taken  to  be 


u{t)  = 


sin(1207rt)  +  sin(3607rt)  +  sin(8007rt) 
0  ,  2/60  <  t  <  10/60 


0  <  <  <  2/60 


with  the  frequencies  again  chosen  so  as  to  strongly  excite  the  first,  second,  fourth  and  sixth 
system  modes  (see  Figure  11).  During  the  final  3/4  of  the  time  interval,  the  voltage  produced 
by  the  patches  cis  well  as  the  acoustic  pressure  at  the  cavity  point  (.6,  .1)  were  calculated  at 
750  uniformly  distributed  times  throughout  the  interval  [5/120,10/60].  These  values  were 
then  used  as  our  simulated  data. 

The  estimated  values  of  the  parameters  parameters  p,  EI,cd1,IC^  and  IC^  are  recorded  in 
Table  4.  In  the  first  four  simulations,  168  cavity  and  19  beam  basis  functions  were  used  to 
calculate  the  data  while  120  cavity  and  15  beam  basis  functions  were  used  in  the  estimation 
of  the  physical  parameters.  The  values  in  columns  3  and  4  were  obtained  with  data  consisting 
of  the  750  voltage  values  generated  by  the  piezoceramic  patches  while  data  consisting  of  both 
voltage  and  acoustic  pressure  values  was  used  to  obtain  the  results  in  columns  5  and  6.  Finally, 
the  results  from  a  simulation  in  which  120  cavity  and  15  beam  baisis  functions  were  used  for 
both  the  generation  of  data  and  estimation  of  parameters  are  reported  in  column  7  of  Table  4. 
The  data  for  this  latter  simulation  consisted  of  the  previously  described  combination  of  voltage 
and  pressure  values. 

It  is  first  noted  that  when  the  same  number  of  basis  functions  are  used  for  the  generation 
of  data  and  estimation  of  parameters,  highly  accurate  results  can  be  obtained  for  a  variety  of 
initial  starting  values  when  using  a  combination  of  voltage  and  pressure  values  as  data  (the 
results  in  column  7  of  Table  4  were  obtained  with  the  initial  values  in  column  5,  Table  i). 
Th.’s  is  consistent  with  the  results  reported  in  Example  2  and  described  in  the  Example  3.  As 
discussed  in  the  latter  example,  the  use  of  a  larger  number  of  beisis  functions  when  determining 
the  data  adds  a  form  of  “numerical  noise”  since  the  data  is  calculated  with  greater  accuracy 
than  are  the  solutions  obtained  during  parameter  estimation.  Hence,  although  no  random 
noise  w^  added  to  the  data,  the  estimated  values  in  columns  3,5  and  6  of  Table  4  differ 
slightly  from  the  “true”  values  listed  in  Table  1  but  are  consistent  with  those  obtained  in  the 
last  example  (see  column  3,  Table  3). 

The  effects  that  the  accuracy  of  the  initial  guesses  had  on  the  estimated  parameters  can 
be  seen  by  comparing  the  results  in  columns  3  and  4  (voltage  data)  as  well  as  in  columns  5 
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and  6  (voltage  and  pressure  data).  The  parameter  values  in  columns  3  and  5  were  obtained 
with  the  highly  accurate  initial  choices  of  column  4,  Table  1  while  those  in  columns  4  and 
6  were  calculated  with  the  less  accurate  values  in  column  5,  Table  1  which  could  reasonably 
have  been  obtained  by  first  visually  fitting  the  data  in  the  frequency  domain.  As  observed,  the 
change  in  starting  values  led  to  fairly  large  changes  in  the  recovered  parameter  values  when 
purely  voltage  data  was  used,  whereas  almost  no  change  was  noted  when  the  data  consisted 
of  voltage  and  pressure  values. 

The  variation  in  the  recovered  parameter  values  obtained  with  voltage  data  can  partially 
be  attributed  to  the  flatness  of  the  criterion  functional  with  strain  observations  as  compared 
to  that  seen  with  acceleration  data.  The  inclusion  of  pressure  values  in  the  data  adds  rich¬ 
ness  as  a  result  of  the  added  information  about  ‘the  acoustic  state  as  well  as  the  fact  that 
higher  frequencies  are  more  easily  observed  in  the  pressure  data  than  in  the  stain  values  (see 
Figure  10).  Due  to  the  added  information  in  the  data  containing  strain  and  pressure  values, 
a  wider  range  of  initial  guesses  could  be  used  since  the  optimization  routine  was  less  like  to 
become  stuck  in  local  minima. 


Siml 

Sim2 

Sim3 

Sim4 

Sim5 

p 

Region  1 

1.3501 

1.3745 

1.3502 

1.3502 

1.3500 

(fc<//m) 

Region  2 

2.1207 

2.1476 

2.1230 

2.1231 

2.1150 

El 

Region  1 

73.178 

76.319 

72.910 

72.908 

73.960 

(Nm^) 

Region  2 

128.022 

123.546 

129.041 

129.051 

1*25.400 

CD  I 

Region  1 

.000992 

.000283 

.000989 

.000988 

.000999 

{Nm^sec) 

Region  2 

.001274 

.000431 

.001295 

.001297 

.001250 

{N 

mfV) 

.007274 

.006709 

.006767 

.006766 

.006700 

IC^iV) 

158.749 

85.832 

172.018 

171.977 

169.999 

Table  4.  Estimated  values  of  p,  El,  cdI,IC^  andK^;  (a)  SimJ:  voltage  data  and  the  initial  guesses 
of  Table  1,  column  5;  (b)  Sim2:  voltage  data  and  the  initial  guesses  of  Table  1,  column  6;  (c)  SimS: 
voltage  and  pressure  data  with  the  initial  guesses  of  Table  1,  column  5;  (d)  Sim^:.  voltage  and  pressure 
data  with  the  initial  guesses  of  Table  1,  column  6;  (e)  SimS:  voltage  and  pressure  data  generate  with 
same  number  of  basis  functions  used  in  parameter  estimation,  initial  guesses  of  Table  1,  column  6. 


Fig.  10.  The  system  response  to  an  initial  periodic  voltage.  The  solid  line  (with  x’s)  denotes  the 
data  containing  10%  noise  while  the  dashed  line  (with  o’s)  is  the  model  response  obtained  with  the 
recovered  parameter  values  in  Table  4,  column  6. 


Example  5:  Patch  Impulse  Excitation,  Patch  Data 

The  input  voltage  in  this  example  weis  taken  to  be  a  narrow  triangle  which  simulates  a 
voltage  impulse  to  the  patches.  As  demonstrated  by  the  strain  and  pressure  plots  in  Figure  12, 
this  causes  the  excitation  of  multiple  system  modes  which  through  time  decay  in  energy  due 
to  the  damping  in  the  beam. 

The  data  was  generated  using  168  cavity  and  19  beam  basis  functions,  and  a  comparison 
was  again  made  between  the  results  obtained  with  data  consisting  solely  of  voltage  values 
(using  the  criterion  functional  in  (4.9))  and  data  made  up  of  both  voltage  and  pressure  mea- 
surenients  (with  the  criterion  functional  in  (4.10)).  In  both  cases,  the  data  wais  calculated  at 
.900  uniformly  spaced  points  throughout  the  time  interval  [1/60, 10/60] 

The  parameter  values  obtained  with  120  cavity  and  15  beam  bcisis  functions  and  strain 
(voltage  values)  data  are  given  in  Table  5.  The  initial  guesses  of  column  5,  Table  1  were 
.used  to  obtain  both  the  results  containing  10%  random  noise  and  those  to  which  no  noise  was 
^ded  (the  noisefree  results  can  be  seen  to  be  quite  close  to  those  obtained  in  the  last  two 


examples).  Finally,  plots  containing  the  time  and  frequency  domain  results  for  the  case  in 
which  10%  random  noise  was  added  to  the  data  are  given  in  Figure  12. 

We  point  out  that  with  this  means  of  system  excitation,  a  much  wider  range  of  initial 
guesses  led  to  convergence  to  the  reported  values  than  was  the  case  in  the  last  example  where 
a  periodic  voltage  was  used  to  excite  the  system.  This  appears  to  be  due  to  the  fact  that  a 
larger  number  of  frequencies  are  excited  and  observed  in  the  strain  data  which  in  turn  leads 
to  “richer”  data.  On  the  other  hand,  the  large  number  of  frequencies  observed  in  the  pressure 
field  (see  Figure  12)  made  it  difficult  to  use  pressure  values  in  combination  with  the  voltage 
as  data  since  an  extremely  good  initial  guess  Wcis  required  in  order  to  obtain  convergence 
of  the  optimization  routine  (in  this  case,  one  must  work  nearly  exclusively  in  the  frequency 
domain  at  first  since  any  deviation  in  the  density  knd  stiffness  parameters  leads  to  frequency 
changes  (see  (5.1))  that  make  estimation  in  the  time  domain  very  difficult).  We  also  noted 
that  similar  problems  were  encountered  when  acceleration  data  was  used  with  this  number  of 
frequencies  excited  due  to  the  fact  that  the  acceleration  provides  a  more  sensitive  measure  of 
beam  movement  than  do  the  patches  which  measure  accumulated  strain.  Hence,  with  a  large 
number  of  frequencies  excited,  the  role  of  the  patches  as  the  sole  sensors  in  the  structural 
acoustic  system  appears  to  improve,  and  in  some  cases  this  method  of  data  observation  may 
be  preferable  since  the  matching  of  a  smaller  number  of  observed  frequencies  may  lead  to  a 
more  tractable  optimization  problem. 


No  Noise 

10%  Noise 

p 

Region  1 

1.3516 

1.3534 

(%/m) 

Region  2 

2.1236 

2.1218 

El 

Region  1 

73.261 

73.969 

Region  2 

128.156 

126.152 

CD  I 

Region  1 

.000993 

.000959 

(Nni^sec) 

Region  2 

.001273 

.001334 

{N 

m/V') 

.006677 

.006188 

K^V) 

173.200 

183.266 

Table  5.  Estimated  values  of  p,  EI,cdI,IC^  and  K.^  using  strain  data;  (a)  No  noise  added  to  the 
data;  (b)  10%  noise  added  to  the  data. 
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Fig.  12.  The  system  response  to  an  initial  voltage  impulse.  The  solid  line  denotes  the  data  containing 
10^  noise  while  the  dashed  line  is  the  model  response  obtained  with  the  recovered  parameter  values 
in  Table  5. 


6  Conclusions  and  Implementation  Issues 

In  this  paper,  a  parameter  estimation  methodology  for  structural  acoustic  systems  in  which 
piezoceramic  patches  on  the  structure  act  as  sensors  and  actuators  has  been  presented.  While 
illustrated  throughout  this  presentation  in  the  context  of  a  2-D  acoustic  cavity  with  a  thin 
beam  at  one  end,  the  methods  developed  here  are  equally  applicable  in  many  3-D  distributed 
parameter  models  (such  as  those  in  [9, 10])  of  actual  experimental  devices  currently  being  used 
to  test  various  modeling,  parameter  estimation  and  control  schemes  involving  piezoceramic 
actuators  and  sensors.  Recently,  these  estimation  methods  have  been  successfully  used  as  an 
integral  part  of  feedback  control  techniques  with  experimental  data  -  these  findings  will  be 
reported  in  a  future  report. 

The  emphasis  in  developing  this  parameter  estimation  methodology  was  on  providing  a 
scheme  which  was  amenable  to  approximation  and  implementation  under  a  variety  of  damping 
eissumptions.  In  doing  so,  conditions  leading  to  well-posedness  and  parameter  convergence 
results  for  both  the  linearized  and  original  nonlinear  problem  were  formulated,  and  suitable 
numerical  techniques  for  approximating  system  dynamics  and  implementing  the  parameter 
estimation  schemes  were  developed. 

To  illustrate  the  method,  several  numerical  examples  illustrating  a  variety  of  techniques  for 
exciting  the  system  and  generating  data  were  presented.  Throughout  the  examples,  emphasis 
was  placed  on  simulating  exterior  forces  and  generating  data  in  a  manner  consistent  with  that 
used  in  actual  structural  acoustics  experiments  and  applications.  This  Wcis  done  so  as  to  gain 
insight  regarding  the  eflFectiveness  of  the  parameter  estimation  method  in  various  settings.  The 
first  step  in  the  examples  was  the  determination  of  the  natural  frequencies  for  the  symmetric 
(in  x)  modes  of  the  coupled  system.  This  knowledge  was  then  used  when  determining  driving 
frequencies  so  as  to  evoke  particular  system  responses  in  the  later  examples.  Moreover,  by 
cornparing  the  results  for  the  system  under  investigation  with  comparable  results  obtained  for 
the  structural  acoustic  system  in  which  no  patches  were  bonded  to  the  beam  as  well  as  analytic 
values  of  the  natural  frequencies  of  the  uncoupled  beam  and  acoustic  wave,  the  qualitative 
effects  of  coupling,  damping,  and  material  changes  due  to  the  presence  of  the  patches  were 
determined. 

A  simulated  acoustic  source  was  used  to  excite  the  beam,  and  hence  the  system,  in  the 
second  example  with  data  consisting  of  acceleration  values  calculated  at  the  center  of  the 
beam.  While  the  numerical  simulations  demonstrated  the  success  of  the  method  for  both 
noisefree  data  and  data  to  which  random  noise  had  been  added,  this  means  of  exciting  the 
system  in  a  manner  that  can  be  accurately  simulated  will  be  difficult  to  implement.  This 
is  due  to  the  fact  that  although  speakers  can  be  used  to  create  multiple  frequency  acoustic 
forces  of  the  type  used  in  the  example,  it  will  be  difficult,  if  not  impossible,  to  cleanly  cut  the 
acoustic  excitation  at  a  given  time  since  the  speaker  and  room  will  continue  to  reverberate 
and  echo  even  after  the  power  is  cut.  This  echo  or  reverberation  will  be  difficult  to  simulate  in 
the  system  model  thus  making  parameter  estimation  through  this  means  of  system  excitation 
difficult  to  use  in  practice. 

A  more  readily  simulated  means  of  exciting  the  system  is  through  the  application  of  a 
prescribed  voltage  to  the  patches  and  this  was  the  input  in  the  remaining  examples  (by  altering 
,  the  frequencies  and  magnitude,  this  voltage  can  be  tailored  so  as  to  evoke  a  desired  system 
response).  In  the  third  and  fourth  examples,  a  numerically  simulated  multiple  frequency. 
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periodic  voltage  was  used  to  initially  excite  the  system  at  which  point  the  voltage  was  cut  and 
the  system  energy  was  allowed  to  begin  decaying  due  to  the  damping  in  the  beam.  In  addition 
to  being  more  accurately  simulated  numerically,  this  also  proved  to  be  a  more  effective  means 
of  exciting  high  frequency  responses  in  the  system  than  was  the  acoustic  excitation.  A  short 
duration  triangular  input  was  applied  to  the  patch  in  the  final  example  to  simulate  an  initial 
voltage  spike. 

In  the  third  example,  acceleration  values  on  the  beam  were  used  as  data,  whereas  voltage 
values  (mecisuring  accumulated  strain  in  the  beam)  and  voltage  values  in  combination  with 
acoustic  pressure  values  were  used  as  data  in  the  fourth  and  fifth  examples.  In  these  examples, 
it  was  found  that  when  only  a  few  system  frequencies  were  excited  (thrc^  or  four),  the  accel¬ 
eration  data  was  sufficiently  sensitive  so  as  to  permit  effective  parameter  estimation  whereas 
highly  accurate  initial  guesses  were  needed  in  order  to  use  beam  data  containing  only  strain 
information  (which  qualitatively  has  the  properties  of  displacement  data).  By  augmenting  the 
strain  data  with  pressure  measurements,  however,  sufficient  information  was  added  so  as  to 
again  lead  to  successful  parameter  identification  with  a  range  of  initial  values.  On  the  other 
hand,  the  excitation  of  large  number  of  system  frequencies  through  the  simulated  voltage 
spike  to  the  patch  led  to  strain  data  which  was  sufficiently  “rich”  so  as  to  permit  success 
of  the  method  while  the  introduction  of  pressure  data  led  to  a  failure  in  the  optimization 
routine  due  to  the  very  large  number  of  frequencies  in  the  pressure  measurements.  Hence  the 
success  of  the  method  with  strain  measurements  comprising  pa^'t  of  the  data  depended  partly 
on  the  number  of  excited  frequencies;  for  a  small  number,  the  strain  measurements  had  to  be 
supplemented  with  acceleration  or  pressure  values  in  order  to  ensure  parameter  convergence 
for  data  containing  any  noise. 
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